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INTRODUCTION 


Octahydro*1 ,3,5,7-tetranitro-l  ,3,5,7-tetrazocine  (HMX)  and  hexahydro- 
1 ,3,5-trinitro-l  ,3,5-s-triazine  (RDX)  (Fig.  1)  have  been  the  focus  of  a  great  deal 
of  research  since  their  discovery  around  the  turn  of  the  century.12  The  acronym 
HMX,  coined  by  researchers  in  Great  Britain,1  derives  from  the  phrase  "High 
Melting  Explosive  of  His  Majesty's  Explosive."  This  code  name  was  based  on 
the  fact  that  in  the  early  syntheses  of  HMX,  RDX  was  present  as  a  major 
contaminant.2  HMX  was  simply  the  higher  melting  of  the  two  explosives.  The 
acronym  RDX,  also  a  code  name  given  by  the  British  researchers,  stands  for 
"Research  Department  Explosive’  or  "Royal  Demolition  Explosive." 
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HMX  RDX 


FIGURE  1 .  Structures  of  HMX  and  RDX 


HMX  and  RDX  are  cyclic  nitramines  which  are  both  high  explosives. 
They  posses  an  explosive  force  approximately  one  and  a  half  times  greater 
than  TNT.1-2  Both  HMX  and  RDX  have  been  used  in  munitions  since  before 
World  War  U,  but  it  has  been  only  in  the  past  two  decades  that  they  have  been 
employed  in  rocket  propellants. 

SOLID  ROCKET  PROPELLANTS 

Solid  rocket  propellants  consist  primarily  of  three  components:  the  fuel, 
the  oxidizer,  and  the  rubbery  binder.3  A  typical  solid  rocket  propellant  is 
composed  principally  of  powdered  aluminum  (Al)  fuel,  ammonium  perchlorate 
(AP)  oxidizer,  and  a  crosslinked  polybutadiene  binder. 

Solid  propellant  performance  is  normally  quoted  in  terms  of  specific 
impulse  (1^).  The  propellant  thrust  (F)  and  the  mass  flow  rate  of  exhaust 
products  (w)  are  related  to  lv  by  equation  1  * 


The  unit  of  F  is  pounds  (force),  and  the  units  of  w  are  pounds  (mass)  per 
second.  The  unit  of  is  normally  given  as  seconds.  However,  the  actual  unit 
is  feet  per  second.  The  l#p  of  a  typical  AI/AP  propellant  ranges  from  260  to  265 
seconds. 

The  effect  of  pressure  on  the  burn  rate  (r)  of  a  propellant  is  given  by 

r  -  a(pc)"  (2) 

where  pc  is  the  pressure  inside  the  combustion  chamber,  a  is  an  empirical 
constant,  and  n  is  the  burning  rate  pressure  exponent.3  The  larger  the  pressure 
exponent  the  greater  the  change  in  propellant  burn  rate  with  increasing 
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pressure.  An  acceptable  value  for  the  pressure  exponent  falls  between  0.2  and 
0.8.  In  developing  a  propellant  formulation  one  attempts  to  maximize  the 
performance  (i.e.  the  1^)  while  maintaining  a  reasonable  pressure  exponent. 

Aluminum/ammonium  perchlorate-based  propellants  produce  copious 
amounts  of  smoke  during  combustion  due  to  the  presence  of  Al203  and  HCI  as 

combustion  products.  This  type  of  highly  visible  exhaust  plume  can  be  a  great 
disadvantage  during  a  combat  situation.  Therefore,  in  recent  years  there  has 
been  a  move  toward  low  signature  or  minimum  smoke  propellants.  Most  of 
these  propellants  incorporate  the  nitramines,  HMX  or  RDX,  as  replacements  for 
both  the  fuel  (Al)  and  the  oxidizer  (AP).  Nitramine-based  propellants  give 
performance  comparable  to  conventional  AI/AP  propellants,  and  their 
combustion  produces  only  gaseous  products  giving  them  almost  invisible 
exhaust  plumes. 

Unfortunately,  nitramine  based  propellants  possess  certain  problems 
which  have  thus  far  limited  their  application.  To  begin  with,  nitramine 
propellants  tend  to  have  large  pressure  exponents,  often  greater  than  0.8.  As  a 
consequence  the  propellant  burn  rate  and  the  combustion  pressure  are  very 
sensitive  to  each  other  such  that  catastrophic  pressure  increases  can  occur  in  a 
few  milliseconds.3 

It  is  difficult  to  vary  the  burn  rate  of  nitramine  propellants  through 
modification  of  the  propellant  formulation.  Techniques  developed  to  modify  the 
burn  rate  of  AI/AP  propellants  such  as  changing  ingredient  particle  size,  adding 
catalysts,  varying  the  amounts  of  ingredients,  etc.,  have  little  effect  on  the  burn 
rate  of  nitramine  propellants.  In  essence  nitramine  propellants  tend  to  burn 
over  a  very  narrow  range  often  referred  to  as  the  "bum  rate  box".  To  make 
matters  worse,  the  burn  rates  of  nitramine  propellants  are  normally  too  slow  for 
most  applications. 
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Al/AP-based  propellants  are  normally  in  the  1.3  hazard  classification 
range.  However,  because  HMX  and  RDX  are  high  explosives,  propellants 
made  from  these  nitramines  tend  to  be  relatively  sensitive,  and  are  normally  in 
the  1.1  hazard  classification  (the  classification  1.1  is  for  compounds  that 
detonate  whereas  the  classification  1 .3  is  for  compounds  that  bum  rapidly  or 
deflagrate).  The  increased  hazard  of  nitramine  propellants  presents  significant 
handling  problems  and  in  general  limits  their  application  to  land  based  systems. 

If  HMX  and  RDX  based  propellants  are  to  be  used  on  a  broad  scale 
significant  improvements  in  the  pressure  exponent,  burn  rate  range,  and 
propellant  hazard  must  be  made.  This,  in  turn,  necessitates  a  fundamental 
understanding  of  the  processes  governing  the  energy  release  in  the 
combustion  of  HMX  and  RDX. 

HMX  AND  RDX  RESEARCH 

The  initial  work  on  HMX  and  RDX  decomposition  was  published  by 
Robertson5  in  1949.  Since  that  time  a  large  body  of  work,  some  of  which  is 
summarized  in  the  reviews,6*15  has  been  devoted  to  studying  all  aspects  of  the 
thermal  processes  (decomposition,  pyrolysis,  deflagration,  combustion,  and 
detonation)  of  HMX  and  RDX.  To  a  large  extent  these  studies  have  simply 
generated  a  myriad  of  proposed  mechanisms,  attesting  to  the  complex  nature  of 
these  thermal  reactions.  Recently,  researchers  have  focused  their  attentions  on 
determining  the  rate  limiting  reactions  in  the  thermal  decomposition  and  high 
pressure  combustion  of  HMX  and  RDX  because  it  is  the  rate  of  the  slowest 
reaction  which  to  a  large  extent  controls  the  global  energy  release  process. 

Kinetic  isotope  effects  provide  a  powerful  technique  for  the  investigation 
of  the  rate  controlling  reactions  of  the  thermal  process  of  energetic  materials 
Recently,  isotope  effects  have  been  used  to  study  the  thermal  decomposition 


and  combustion  of  HMX  and  RDX.16*20  The  results  of  these  experiments 
implicate  a  C*H  bond  rupture  as  the  global  rate  controlling  phenomena  in  these 
thermochemical  reactions. 

The  isotope  effects  studies  have  important  implications  in  the 
improvement  of  nitramine-based  propellants.  However,  the  relevance  of  results 
for  pure  HMX  and  RDX  is  limited  because  of  the  possibly  significant  influence 
that  energetic  propellant  ingredients  may  exert  on  the  mechanism  of  nitramine 
propellant  combustion. 

In  a  continuation  of  previous  combustion  work  on  HMX16  and  RDX17, 1 
will  report  isotope  effects  in  the  combustion  of  model  nitramine  propellants. 
These  results  will  be  compared  to  the  results  for  the  combustion  of  the  pure 
nitramines  in  a  effort  to  draw  some  conclusions  about  the  effect  of  propellant 
ingredients  on  the  nitramine  combustion  mechanism.  Finally,  possible 
explanations  will  be  offered  for  the  observed  isotope  effects  and  their  kinetic 
implications. 
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KINETIC  DEUTERIUM  ISOTOPE  EFFECT  THEORY 


Kinetic  isotope  effects  are  the  primary  tool  used  in  this  research  to  study 
the  combustion  kinetics  of  nitramines  and  nitramine  based  propellants.  Thus,  it 
is  important  to  establish  a  clear  theoretical  understanding  of  kinetic  isotope 
effects  before  conclusions  can  be  made  from  our  data.  Furthermore,  a 
theoretical  foundation  is  needed  such  that  a  critical  evaluation  of  previous 
isotcpe  effect  studies  concerning  the  decomposition  and  combustion  kinetics  of 
energetic  materials  can  be  made. 

A  kinetic  isotope  effect  refers  to  the  change  in  the  rate  at  which  a 
molecule  undergoes  a  specific  reaction  when  an  atom  (or  atoms)  in  the 
molecule  is  (are)  replaced  with  a  different  isotope.  Isotope  effects  are  normally 
expressed  as  a  ratio  of  rate  constants.  In  the  case  of  the  substitution  of 
deuterium  for  protium,  the  isotope  effect  is  most  often  represented  as  l^/ko, 
where  kH  is  the  rate  constant  for  the  molecule  containing  protium  and  kQ  is  the 
rate  constant  for  the  molecule  containing  deuterium. 

According  to  the  Born-Oppenheimer  approximation,  the  potential  energy 
surface  of  a  molecule  depends  only  on  the  attractions  and  repulsions  of 
electronic  and  nuclear  charges,  not  on  nuclear  mass.2122  Nuclear  motion,  on 
the  other  hand,  does  depend  on  the  mass  of  the  nucleus,  and  in  turn  the  relative 
nuclear  motion  determines  the  vibrational  energy  content  of  the  molecule.23 
Consequently,  molecules  differing  only  by  isotopic  substitution  move  on  the 
same  potential  energy  surface  but  differ  in  vibrational  energy  content.  To  a 
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large  extent,  it  is  this  vibrational  energy  difference  which  gives  rise  to  kinetic 
isotope  effects. 

Kinetic  isotope  effects  have  long  been  recognized  for  their  potential  in 
the  study  of  reaction  rates  and  mechanisms.  Normal  substituent  effects  do  not 
always  give  a  clear  understanding  of  a  system  due  to  changes  in  the  potential 
energy  surface  and  possible  mechanistic  changes  which  occur  when  going 
from  one  substituent  to  another.  The  potential  energy  surface  and  the  reaction 
mechanisms  of  molecules  differing  only  by  isotopic  substitution  are  essentially 
the  same.  This  allows  kinetic  isotope  effects  to  be  used  to  investigate  the  nature 
of  a  single  potential  energy  surface  and  thus  a  single  transition  state.23 

Isotope  effects  are  categorized  as  being  either  primary  or  secondary.  A 
primary  (1°)  isotope  effect  occurs  when  a  bond  to  the  isotopic  atom  is  broken 
during  the  rate  limiting  step.  Secondary  (2°)  isotope  effects  are  those  effects 
which  do  not  involve  the  rupture  of  a  bond  to  the  isotopic  atom.  Primary  isotope 
effects  are  generally  larger  than  secondary  effects.  However,  as  we  shall  see, 
determining  whether  an  observed  isotope  effect  is  primary  or  secondary  is  far 
from  trivial. 

For  single  step  reactions,  kinetic  isotope  effects  are  observed  when 
bonds  to  an  isotopic  atom  undergo  changes  when  going  to  the  activated 
complex.  In  the  case  of  multi-step  reactions,  isotope  effects  will  be  seen  if 
bonds  to  an  isotopic  atom  undergo  changes  In  a  slow  step  of  the  reaction 
mechanism. 

ORIGINS  OF  KINETIC  ISOTOPE  EFFECTS 

The  theoretical  foundation  for  the  treatment  of  kinetic  isotope  effects  was 
first  advanced  by  Bigeleisen  and  Mayer  In  1947.24  Kinetic  isotope  effect  theory 
is  based  on  the  "transition  state  theory*  of  Eyring25  which  assumes  that 


reactants  are  in  thermal  equilibrium  with  an  "activated  complex"  or  "transition 
state."  Furthermore,  it  is  this  activated  complex  which  decomposes  in  the  rate 
limiting  step  of  the  overall  process  into  products.26  Bigeleisen  and  Mayer 
employed  the  tenets  of  transition  state  theory  coupled  with  statistical  mechanics 
to  derive  relationships  which  provide  a  reasonable  explanation  for  observed 
kinetic  Isotope  effects 

The  following  discussion  focuses  on  the  isotope  effects  arising  when 
deuterium  is  substituted  for  profium,  hereafter  referred  to  as  Kinetic  Deuterium 
Isotope  Effects  (KDIE).  For  more  detailed  information  the  reader  is  directed  to 
the  reviews23'26-32  listed  in  the  reference  section  of  this  thesis 


Zero  Point  Energy 

A  great  majority  of  kinetic  isotope  effects  result  from  differences  in  zero 
point  energy  between  normal  and  isotopically  labeled  molecules. 
Consequently,  a  brief  explanation  of  the  origin  of  zero  point  energy  differences 
is  necessary  for  a  complete  understanding  of  kinetic  isotope  effects. 

For  a  molecule  A-H  where  H  is  hydrogen  and  A  is  the  remainder  of  the 
molecule,  the  energy  levels  for  the  vibrations  of  A-H  will  be  given  by 


E 


hv 


(3) 


where  h  is  Planck's  constant,  n  is  the  vibrational  quantum,  and  v  is  the  A-H 


vibrational  frequency.  This  of  course  assumes  that  the  A-H  stretch  is  an 
isolated  vibration  and  not  part  of  a  normal  mode.  The  lowest  energy  level  (n  « 
0)  is  the  zero  point  energy  (Eq);3334  it  is  given  by 


E0 


(4) 
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At  normal  temperatures  almost  all  molecules  occupy  the  n  -  0  vibrational 
energy  level.33-35 

Assuming  A-H  behaves  like  a  harmonic  oscillator,  the  vibrational 
frequency  will  be  given  by 


where  k  is  the  force  constant  for  the  particular  vibrational  motion,  and  p  is  the 
reduced  mass  of  A-H.  If  A  is  significantly  larger  than  H,  the  infinite  mass 
approximation22  can  be  used  to  obtain 


where  mH  is  the  mass  of  the  hydrogen  atom. 

Combining  equations  4  and  6  we  obtain  the  following  expression  for  the 
zero  point  energy  of  a  vibration. 

6° -£75  ,7> 

Molecules  differing  only  in  isotopic  mass,  such  as  A-H  versus  A-D 
where  D  is  deuterium,  move  along  the  same  potential  energy  surface.  Since 
the  vibrational  force  constants  depend  solely  on  this  potential  energy  surface, 
they  are  the  same  for  both  A-H  and  A-D.  Thus,  with  this  in  mind  equation  7  can 
be  reduced  to 

E°-7k  (8) 

where  m  is  the  mass  of  the  light  atom  (i.e.,  H  or  D).  Consequently,  the 
vibrations  of  A-D  will  have  lower  zero  point  energies  than  the  vibrations  of  A-H. 

KDIE  in  a  Hydrogen  Atom  Transfer 

In  an  effort  to  make  a  complex  subject  more  understandable,  further 
discussions  of  kinetic  isotope  theory  will  center  on  the  isotope  effects  arising  in 


a  simp  to  hydrogen  transfer  reaction.  The  logic  used  in  the  following  discussion 
wiH  be  generally  applicable  to  most  primary  and  secondary  KDIEs. 

Consider  the  following  reaction  involving  the  transfer  of  a  pro ti urn  atom 
(H),  or  a  deuterium  atom  (0)  from  molecule  A  to  molecule  B.26 

AH(D)  ♦  B  ^  A-H(D)-B  — ►  A  ♦  BH(D)  (9) 

(The  discussion  below  will  focus  primarily  on  the  reaction  involving  protium  with 
the  understanding  that  what  is  said  is  also  true  for  deuterium.)  Following 
transition  state  theory,  AH  and  B  are  in  thermal  equilibrium  with  the  activated 
complex  A*"*H***‘B  which  subsequently  decomposes  into  the  products  A  and 
BH.  Figure  2a26-31  32  shows  a  representation  of  a  potential  energy  surface  for 
this  reaction.  The  BH  bond  length  is  plotted  versus  the  AH  bond  length. 
Energy  is  plotted  perpendicular  to  the  plane  of  the  paper  with  the  contours 
connecting  points  of  equal  energy.  The  outermost  contours  are  the  highest  in 
energy  with  the  energy  of  the  contours  decreasing  towards  the  center.  The  line 
of  minimum  energy  along  which  the  transition  state  lies  is  represented  by  a 
dashed  fine  and  is  the  reaction  path. 

Figure  2b31-32  represents  the  reaction  path  from  Figure  2a  with  cross 
sections  of  the  potential  energy  surface  corresponding  to  the  reactants, 
activated  complex,  and  products.  All  three  potential  wells  show  positive 
curvature  and  thus  positive  force  constants  perpendicular  to  the  reaction  path. 
In  the  case  of  the  activated  complex,  the  potential  well  corresponds  to  a 
symmetric  stretch  (i.e.,  «-A’"’H‘"*B-»). 

The  potential  energy  surface  (Fig.  2b)  shows  a  neutral  curvature  in  the 
direction  of  the  reaction  path  for  both  the  reactants  and  products,  indicating  a 
zero  force  constant  for  motion  in  that  direction.  Thus,  there  is  no  resistance  to 
the  extension  of  the  AH-B  distance  for  the  reactants  and  the  A-HB  distance 


10 


FIGURE  2a?8,31,32  Potential  Energy 
Surface  for  the  reaction  in 
equation  7.  The  dashed  line 
is  the  reaction  path. 


(3) 


A  +  BH(D) 


FIGURE  2b31 32  Reaction  Path 
from  Figure  2a.  Potential  Energy 
Surface  cross  sections  are  from 
the  locations  shown  in  Figure  1  a 
and  are  for  the  following:  (1) 
reactants;  (2)  activated  complex; 
(3)  products. 


A  •••  BH(D) 


for  the  products.  At  the  activated  complex  the  potential  energy  surface  shows  a 
negative  curvature  along  the  reaction  path.  Consequently,  motion  of  the 
activated  complex  along  the  reaction  path,  which  corresponds  to  an  asymmetric 
stretch  (i.e.,  «-A***,H-»«-B  and  A-»«-H"“B-»),  is  unstable.  This  follows 
directly  from  transition  state  theory,  which  requlm  ^rational  mode  to 

become  translation  along  th^  reaction  path  26'3® 

The  total  zero  point  energy  of  a  moleo*  m  the  sum  of  the 

zero  point  energies  (j  hv)  of  all  the  vibrational  i  esc  a*  alecule.  In  the 

case  of  the  product  (A  and  BH)  and  the  reactai  '  ^  >  molecules,  the 

total  zero  point  energies  are  the  result  of  the  oumnoutions  of  the  3N-6 

vibrational  modes  of  these  molecules,  where  N  is  the  number  of  atoms  in  the 
molecule.  The  total  zero  point  energy  of  the  activated  complex  (A****H*— B)  is 
made  up  of  the  contributions  of  3N-7  vibrational  modes  (one  mode  is  lost  to 
decomposition  along  the  reaction  path).21*26  The  solid  line  across  the  bottom 
portions  of  the  three  potential  wells  in  Figure  1b  represents  the  zero  point 
energy  level  for  the  protiated  molecules,  and  the  dashed  line  represents  the 
zero  point  energy  level  for  the  deuteriated  molecules. 

According  to  transition  state  theory,  the  rate  constant  k  for  a  reaction  is 
given  by 

KJ  - 

k--|-KK  (10) 

h 

where  kB  is  the  Boltzmann  constant,  T  is  temperature,  and  h  is  Planck's 
constant.  The  transmission  coefficient  (k)  represents  the  probability  that  a  given 
system  moving  in  the  direction  of  the  reaction  path  will  in  fact  pass  over  the 
potential  energy  barrier  to  the  products  side.23  The  quantity  K  is  the  equilibrium 
constant  between  the  reactants  and  activated  complex.  For  the  hydrogen 
transfer  reaction  it  takes  the  form  of 
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(11) 
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AH^B 


exp[-|r] 


where  Qj^B  is  the  partition  function  for  the  activated  complex,  QAH  and  QB  are 

the  molecular  partition  functions  for  the  reactants,  and  R  is  the  universal  gas 
constant  (the  Cross  of  Lorraine  ($)  is  used  to  identify  quantities  associated  with 
the  activated  complex/transition  state).  The  activation  energy  for  the  reaction 
(£,)  is  simply  the  difference  between  the  zero  point  energy  of  the  reactants  (E0) 
and  the  zero  point  energy  of  the  activated  complex  (E£). 

Combining  equations  10  and  11,  the  following  expression  for  the  rate 
constant  of  the  hydrogen  transfer  reaction  is  obtained. 


IC  Vir  QaHB  nirnf  E* 


(12) 


As  was  pointed  out  earlier,  KOIEs  are  represented  as  ratios  of  rate 
constants  (i.e.,  kn/ko)-  Thus,  if  we  take  the  ratio  of  the  rate  constant  expressions 
for  the  protium  and  deuterium  reactions  we  generate  the  following  expression 
for  the  KDIE. 
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This  simplifies  to  give 


kH  _  kh  ^ahb^ad 

*0  *D  Q*  Q 
WADB  AH 


l  RT  J 


(14) 


The  quantity  [AE  is  merely  the  difference  between  the  activation  energy  for 

the  reaction  with  protium  (Ea  h)  and  the  activation  energy  for  th9  reaction  with 
deuterium  (Ea  D). 
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Thera  are  essentially  two  parts  to  equation  14,  each  of  which  may 
contribute  to  the  magnitude  of  the  KDIE.  The  first  part  is  the  preexponential;  it 
consists  of  the  ratio  of  the  transmission  coefficients  and  the  ratio  of  the  partition 
functions.  Under  most  circumstances  both  of  these  ratios  are  very  close  to  unity. 
As  such  they  are  normally  not  considered  when  evaluating  KDIEs.  However, 
there  are  some  instances,  which  will  be  discussed  later,  where  the  observed 
isotope  effects  are  due  almost  entirely  to  the  preexponential. 

The  second  portion  of  equation  14  ftAE,]£)  is  the  zero  point  energy 
|  contribution.  The  value  of  (AEJ$,  in  general,  determines  the  magnitude  of  most 

observed  isotope  effects.  In  turn,  the  size  of  [AEa]$  depends  solely  on  the 
changes  in  zero  point  energy  which  occur  when  going  from  reactants  to 
activated  complex. 

To  understand  how  [AEJ$  arises  we  will  look  back  at  the  hydrogen 
transfer  reaction.  As  stated  earlier,  the  total  zero  point  energy  for  a  molecule  is 
made  up  of  the  contributions  from  all  the  vibrational  modes  of  the  molecule. 
More  importantly,  the  zero  point  energies  of  the  vibrational  modes  of  the 
protiated  and  deuteriated  molecules  are  the  same  except  for  those  vibrations 
involving  the  isotopic  atom.  Therefore,  the  difference  in  the  total  zero  point 
energies  between  the  protiated  and  deuteriated  molecules  in  the  reactants 
((AE<j{J)  and  activated  complex  ((AE£)$)  is  due  only  to  those  vibrations  involving 

an  isotopic  atom. 

From  Figure  3,  which  is  a  more  detailed  representation  of  the  potential 
curves  for  the  reactants  and  activated  complex  from  Figure  2b,  it  should  be  fairly 
clear  that  if  the  zero  point  energy  difference  in  the  reactants  is  the  same  as  the 
zero  point  energy  difference  in  the  activated  complex  (i.e.,  (AE0]£  *  IaE^Jq  ), 

then  the  activation  energies  for  the  transfer  of  protium  and  the  transfer  of 
deuterium  should  be  the  same  (i.e.,  Ea  H  ■  Ea  0).  H  this  were  the  case,  then 
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[AEJq  would  be  zero,  and  consequently  there  would  be  no  zero  point  energy 

contribution  to  the  KDIE.  If,  however,  the  zero  point  energy  differences  are  not 
the  same  (i.e.,  [AE0]g  #  (AEfcjg  )• thd  activation  energies  will  differ  (i.e.  Ea  H  * 
E^d),  resulting  in  a  non-zero  value  for  [AE,]|j  and  a  net  kinetic  isotope  effect. 


FIGURE  3? 1,36  Potential  curves  for  the 
the  reactants  and  activated  complex. 

For  [AEJq  to  have  a  non-zero  value,  the  zero  point  energy  of  at  least  one 
vibration  involving  an  isotopic  atom  must  change  when  going  from  reactants  to 
activated  complex.  One  can  see  from  equation  5  that  a  change  in  zero  point 
energy  (E0)  requires  a  change  in  the  vibrational  force  constant  (k). 

To  illustrate  the  effect  of  force  constant  changes,  I  will  examine  a  KDIE  in 
a  hydrogen  transfer  reaction  (eq.  9)  which  arises  from  a  change  in  force  constant 
when  the  stretching  vibration  of  A-H  in  the  ground  state  becomes  the  symmetric 
stretch  in  the  activated  complex  (i.e.,  «-A’"*H— Suppose  the  force 
constant  for  the  A-H  stretch  decreases  from  465  N/m  for  the  ground  state  to  1 16 
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N/m  tor  the  activated  complex.  This  corresponds  to  lowering  the  A*H  frequency 
from  2800  cm'1  to  1400  cm*'  (the  A-D  frequency  changes  form  1979  cm'1  to  989 
cm*1)  for  the  ground  state  and  activated  complex,  respectively  (the  value  for  the 
activated  complex  is  the  frequency  of  the  hydrogen  in  the  symmetric  stretch).  For 
this  change  in  force  constant,  the  zero  point  energy  for  the  A-H  stretch 
decreases  from  E0  H  -  16.74  kJ/mol  to  E|>H  -  8.37  kJ/mol,  and  the  zero  point 
energy  for  the  A*D  stretch  decreases  from  E0  0  -  11  80  kJ/mol  to  E|  0  -  5.90 

kJ/mol.  This,  in  turn,  results  in  a  lowering  of  the  zero  point  energy  differences 
from  (AEolg  «  4.94  kJ/mol  for  the  reactants  to  [AEfclg  -  2.47  kJ/mol  for  the 
activated  complex.  The  difference  in  these  two  values  (note:  (aE01d  -  [aE$j  - 
[AEJB)  gives  [AEjg  -  -2.47  kJ/mol.  If  you  place  this  value  into  equation  1 4  and 
assume  the  preexponential  to  be  unity,  then  the  resulting  KDIE  would  be  2.7  at 
298  K. 

From  the  above  example  one  can  see  that  a  decrease  in  the  force 
constant  when  going  from  reactants  to  activated  complex  will  give  a  value  of 
kn/kp  greater  than  one.  Conversely,  an  increase  in  the  force  constant  will  result 
in  a  value  of  kH/ko  less  than  one. 

It  must  pointed  out  that  the  above  example  dealt  witn  only  one  hydrogen 
atom,  one  vibration,  and  one  force  constant.  In  real  systems  there  may  be 
numerous  vibrational  modes  involving  isotopic  atoms,  which  undergo  changes 
in  their  force  constants.  Each  one  of  these  vibrations,  in  turn,  would  contribute  to 
the  overall  magnitude  of  [AEJq  and  thus  contribute  to  the  value  of  k^/kQ. 

TYPES  OF  PRIMARY  KDIE 

KDIEs  can  be  categorized  into  different  types  according  to  the 
characteristics  of  their  activated  complex/transition  state.  One  of  the  first 
attempts  to  classify  primary  effects  was  made  by  Westheimer26  who 


16 


distinguished  two  different  types  of  linear  transition  states.  More  recently 
Kwart,2*  incorporating  Westheimer's  categories,  classified  primary  KDIE  into  four 
different  types.  Each  of  these  types  will  be  discussed  in  detail  below. 

To  diverge  momentarily,  equation  14  is  fine  for  discussing  tne  theoretical 
aspects  of  KDIE.  However,  it  is  of  limited  use  when  explaining  experimental 
results.  Kinetic  data  are  normally  obtained  from  experiments  performed  at 
several  temperatures.  These  data  are  then  fitted  to  the  Arrhenius  relation 
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to  obtain  the  activation  energy  (Ea)  and  preexponential  or  frequency  factor  (A). 
The  Arrhenius  equation  is  an  empirical  relationship;  however,  there  is  some 
theoretical  significance  to  both  E,  and  A.37 

If  one  performs  an  operation  on  equation  1 5  similar  to  what  was  done  to 
generate  equation  14,  one  obtains  the  following  expression  for  the  KDIE.36 
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This  equation  is  very  similar  to  the  relationship  for  the  KDIE  derived  Irom 
transition  state  theory  (eq.  14).  In  fact,  for  this  level  of  qualitative  discussion  it  is 
reasonable  to  say  that  the  ratio  of  the  Arrhenius  preexponentials  is 
approximately  equal  to  the  preexponential  portion  of  equation  14. 
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Because  equation  16  is  generated  from  experimental  results,  it  will  be  used  in 
the  following  discussion  of  the  four  general  types  of  KDIE. 

Three  out  of  the  four  categories  of  KDIE  involve  a  linear  transition  state  of 
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the  type  shown  in  equation  9  for  the  transfer  of  a  hydrogen  atom  from  A  to  B. 
The  first  category  involves  a  linear  transition  state  where  the  stretching  vibration 
is  truly  symmetrical.26  A  and  B  move  in  opposite  directions 

while  H  is  motionless.  Consequently,  the  frequency  and  the  zero  point  energy  of 
this  symmetric  stretch  are  independent  of  whether  H  or  D  is  present  (i.e.,  [AE^|J 

>  0).  A  transition  state  of  this  type  gives  a  relatively  large  KDIE  because  the 
activation  energy  difference  is  equal  to  the  difference  in  zero  point  energies  for 
the  reactants  ( [AEj|J « (AE0]|J  ).36  Assuming  a  stretching  frequency  for  A-H  of 
3000  cm*1,  the  value  of  [AEjjJ  would  be  -4.65  kJ/mol,  and  the  observed  KDIE 
would  be  around  6.5  at  298  K.  The  values  of  A^/Aq  for  the  symmetrical  activated 
complex  extend  over  a  narrow  range  from  0.7  to  1 .2.38 

The  second  category  of  KDIE  involves  an  unsymmetrical  linear  transition 
state.36  This  means  that  H,  in  the  activated  complex,  is  bonded  more  strongly  to 
either  A  orB.  Thus,  when  A  and  B  movj  in  opposite  directions  in  the 
"symmetric"  stretch  («-A“"H""B->),  H  will  follow  the  molecule  it  is  bonded  to 
more  strongly.  Consequently,  the  frequency  and  in  turn  the  zero  point  energy  of 
this  so  called  "symmetric"  stretch  will  depend  on  whether  H  or  D  is  present.  As 
can  be  seen  from  the  example  given  in  the  latter  part  of  the  previous  section, 
when  there  is  motion  of  the  isotopic  atom  in  the  activated  complex  [AE^]|J  *■  0. 
Because  [AE0)q  -  [AE  ^  «  [AEa]£,  the  activation  energy  differences  and 
consequently  the  l^/kg  values  will  be  less  than  for  the  truly  symmetrical  case.  In 
practice,  KDIEs  attributed  to  unsymmetrical  linear  transition  states  have  values 
between  2  and  5  at  298  K  and  have  values  of  Ah/Ad  between  0.7  and  1 .2.36-38 

A  third  type  of  KDIE,  possessing  a  linear  transition  state,  depends  on  the 
wave  characteristics  of  the  hydrogen  atom.  The  Heisenberg  uncertainty 
principle  allows  delocalization  of  a  particle  across  a  space  corresponding  to  its 
deBroglie  wavelength.21-39  If  the  potential  barrier  is  narrow  enough  that  the 
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distance  between  A  and  B  across  the  barrier  is  equal  to  the  deBroglie 
wavelength  of  hydrogen,  then  the  hydrogen  atom  can  travel  across  the  barrier  to 
the  products  side  even  if  its  energy  is  below  that  of  the  activated  complex.  This 
process  is  referred  to  as  quantum  mechanical  tunneling,  and  it  produces  kinetic 
rate  constants  which  are  much  larger  than  what  is  expected  by  normal  transition 
state  theory.21'32  Because  the  wavelength  of  protium  is  longer  than  that  for 
deuterium,  departures  from  classical  behavior  will  be  greatest  for  protium. 
Consequently,  KDIEs  for  reactions  involving  tunneling  are  normally  quite  large; 
Kh/Kq  is  often  greater  than  20.32  When  applying  equation  16,  which  is  based  on 
classical  transition  state  theory,  to  the  nonclassical  process  of  tunneling,  the 
values  of  [AEJJJ  will  appear  to  be  much  less  than  [AE|]g,  which  in  turn  results  in 

values  of  Ah/Aq  «  O.7.38 

An  alternative  to  the  three  types  of  linear  transition  states  discussed 
above  is  a  non-linear  hydrogen  atom  transfer,  an  example  of  which  is  shown 
below. 


The  non-linear  or  bent  transition  state  allows  relatively  free  movement  of  the 
hydrogen  atom.  This  is  in  contrast  to  a  linear  transition  state  where  the 
stretching  vibrations  of  the  hydrogen  atom  are  limited  by  the  relatively  close 
proximity  of  molecules  A  and  B.36  Stretching  vibrations  involving  H  change  little 
when  moving  from  reactants  to  a  non-linear  transition  state.  Therefore,  the 
stretching  vibrations  do  not  contribute  significantly  to  the  value  of  [AEJ$. 

In  a  non-linear  transition  state  the  bending  modes  of  H  are  undergoing 
the  most  change.  However,  since  the  energies  of  the  bending  modes  are 
significantly  less  than  for  stretching  modes,  they  too  contribute  very  little  to  the 
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value  of  {AEjjj.36  Overall,  for  a  bent  transition  state  (AE|]g  is  almost  the  same  as 
(AEdjj,  thus  giving  a  small  value  for  the  difference  in  activation  energies  ([aEJq) 
Because  of  the  limited  contribution  from  zero  point  energy,  KDIEs  for  reactions 
with  non-linear  transition  states  are  due  almost  entirely  to  the  preexponential 
ratio  Ah/Aq.  Non-linear  transition  states  have  KDIE  values  around  2  at  298  K 
and  can  be  distinguished  from  unsymmetricai  linear  transition  states  by  values  of 
A„^A0>1.4.3C 

SECONDARY  KDIE 

Secondary  KDIEs  encompass  all  those  isotope  effects  which  do  not 
involve  the  breaking  of  a  bond  to  the  isotopic  atom  in  the  rate  limiting  step  of  the 
reaction.  Secondary  KDIEs,  like  primary,  arise  from  changes  in  vibrational  force 
constants  of  the  normal  modes  involving  the  isotopically  substituted  bonds 
between  the  reactants  and  transition  state.40  In  other  words,  secondary  KDIEs 
are  due  solely  to  zero  point  energy  contributions  (i.e.,  [AE  J{J).41-42 

Halevi43  has  classified  secondary  KDIEs  into  two  different  types  according 
to  whether  or  not  there  Is  a  spatial  reorientation  of  the  bond  to  the  isotopic  atom 
when  going  from  reactants  to  activated  complex.  Secondary  isotope  effects  of 
the  first  kind  involve  spatial  reorientation,  whereas  secondary  isotopes  of  the 
second  kind  do  not. 

Secondary  isotope  effects  of  the  first  kind  have  generally  been  attributed 
to  the  rehybridization  of  the  atom  the  isotope  is  bonded  to.27  43  This  requires  that 
the  isotopic  atom  be  alpha  to  the  reaction  coordinate  (i.e.,  alpha  to  the  bond 
breaking). 

Streitwieser  et  al.44  modeled  the  secondary  KDIE  that  would  occur  when 
a  carbon  atom  rehybridizes  from  sp3  to  sp2  as  shown  in  the  following  reaction 
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To  calculate  the  secondary  isotope  effect,  Streitwieser  et  at.  used  the  following 
expression  derived  from  transition  state  theory. 

K  t*t  f  0.1865  (v  _v*T 

k~  •  1  lexpL — t — l  Hi  HiJJ  (20) 

D  j 

To  arrive  at  this  simplified  expression,  only  zero  point  energy  contributions  were 
taken  into  account,  and  the  approximation  vD  »  v^l  .35  was  used. 

In  calculating  the  secondary  isotope  effect,  Streitwieser  et  al.44  looked  at 
changes  in  the  zero  point  energies  of  only  three  vibrational  frequencies,  one 
stretching  and  two  bending.  To  calculate  KDIE  from  equation  20  requires  values 
for  vibrational  frequencies  in  both  the  ground  state  (vHi)  and  the  activated 
complex  (vjlj,).  Values  for  the  C-H  stretching  and  bending  frequencies  of  the 
ground  state  sp3  carbon  were  easily  assigned.  However,  obtaining  the 
frequencies  for  the  three  vibrations  in  the  activated  complex  required  much  more 
effort.  Using  several  different  C-H  vibrations  of  sp2  carbon  atoms  as  models, 
Streitwieser  et  al.  made  a  best  estimate  of  the  values  of  v^. 

Figure  422-44  shows  a  breakdown  of  the  contributions  of  each  of  the  three 
vibrational  modes  to  the  KDIE.  The  greatest  change  in  frequency  occurs  when 
going  to  the  out-of-plane  bending  vibration  in  the  transition  state,  and  as  such 
this  vibration  makes  the  largest  contribution  to  the  KDIE.  Since  equation  18 
requires  that  the  product  of  each  vibrational  contribution  be  taken,  the  resulting 
value  for  the  KDIE  from  the  data  in  Figure  4  is  1 .48  at  298  K. 
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FIGURE  A.2  24  Calculation  of  Secondary  Isotope 
Effects  at  298  K  for  Rehybridization  of  Carbon 
from  sp3  to  sp2. 


Experimentally,  values  for  secondary  effects  of  the  first  kind  are  in  the 
range  of  1.15  to  1.25  per  deuterium  at  room  temperature.22  For  a  reaction  of  the 
type  in  equation  19,  Streitwieser  observed  a  KDIE  of  only  1.15.44  The  difference 
in  the  experimental  (1.15)  versus  the  calculated  (1.48)  KDIE  values  above  may 
indicate  that  the  hybridization  in  the  activated  complex  lies  somewhere  between 
sp3  and  sp2.  This  would  result  in  a  smaller  change  in  the  C-H  bending 
frequencies  when  going  to  the  transition  state  and  thus  give  a  smaller  KDIE. 

Secondary  isotope  effects  of  the  second  kind  arise  when  the  isotopic 
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atom  is  greater  than  one  bond  away  from  the  reaction  coordinate  (i.e.,  (5,  y,  8). 
Although  there  is  no  spatial  reorientation  of  bonds  to  the  isotopic  atom  in  this 
type  of  secondary  effect,  there  are  still  changes  in  the  force  constants  of  these 
bonds  when  going  to  the  activated  complex.  The  cause  of  secondary  isotope 
effects  of  the  second  kind  has  been  much  debated.23  43  44  They  have  been 
attributed  to  hyperconjugation,  inductive  effects,  and  steric  effects.  It  seems  in 
most  instances  that  hyperconjugation,  which  is  the  interaction  of  the  empty  p- 
orbitals  of  the  carbonium  ion  with  the  adjacent  carbon-hydrogen  o  bond,  is  the 
primary  cause  for  this  type  of  secondary  effect.  This  is  supported  by  the  fact  that 
in  cases  of  multiple  isotopic  substitution,  such  as  a  P  CD3,  there  appears  to  be  a 
conformation  dependence  of  the  size  of  the  secondary  KDIE.  There  are 
instances,  however,  where  both  steric  and  inductive  effects  may  play  an 
important  role.  The  size  of  secondary  isotope  effects  of  the  second  kind  is 
normally  in  the  range  of  1 .1  to  1 .2  per  deuterium  at  298  K.23 

To  conclude,  secondary  KDlEs  are  small,  primarily  because  frequency 
changes  are  small.  Bigeleisen  and  Wolfsberg22  calculated  a  secondary  isotope 
effect  considering  only  changes  in  stretching  frequencies.  For  a  decrease  in  the 
stretching  force  constant  by  a  factor  of  two,  they  calculated  a  secondary  KDIE  for 
a  single  deuterium  of  1 .74.  Secondary  KDIE  values  this  large  are  never  realized 
experimentally  because  the  vibrations  undergoing  the  most  change  in  their  force 
constants  are  bending  and  not  stretching  vibrations.  Normally,  experimental 
values  for  secondary  KDIEs  at  room  temperature  will  range  from  1  0  to  1 .25. 

TEMPERATURE  DEPENDENCE  OF  KDIE 

The  temperature  dependence  of  KDIE  may  be  understood  qualitatively  by 
looking  back  at  equation  16.  The  preexponential  ratio  (AH/AD)  can  be 
considered  to  be  relatively  temperature  independent,36  whereas  the  exponential 
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HAEJ$/RT)  is  quite  obviously  temperature  dependent.  We  know  from  the 
above  discussion  that  large  KDIE  (with  the  exception  of  tunneling)  are  due 
almost  entirely  to  zero  point  energy  contributions.  As  such  KDIE  of  this  type 
should  exhibit  a  profound  dependence  on  temperature.  In  contrast,  KDIE  with 
only  minimal  zero  point  energy  contribution  should  change  little  with 
temperature. 

The  maximum  values  of  [AEJ$  occur  for  reactions  with  linear  symmetrical 
transition  states.23’26  27-32  Consequently,  KDIEs  for  this  type  of  reaction  will  vary 
the  most  with  temperature.  For  example,  using  equation  16  and  assuming  only 
zero  point  energy  contributions,  a  KDIE  of  7.0  at  298  K  would  become  a  KDIE  of 
3.4  at  473  K. 

For  linear  unsymmetrical  transition  states,  the  values  of  [AE  are  smaller 
and  consequently  the  temperature  dependences  of  the  KDIEs  are  less  than  for 
the  symmetrical  case.  Using  equation  16  in  the  same  manner  as  above,  a  KDIE 
for  the  unsymmetrical  case  that  is  2.0  at  298  K  would  decrease  only  to  1 .55  at 
473  K. 

Non-linear  transition  states  have  values  of  [AEJ$  that  are  very  close  to 
zero.  As  such,  the  preexponential  AH/A0is  the  primary  source  of  the  KDIE  for 
reactions  of  this  type.  Because  An/AD  varies  little  with  temperature,  KDIEs  for 
reactions  with  non-linear  transition  states  tend  to  be  temperature  independent. 
In  fact  a  constant  KDIE  over  temperature  ranges  greater  than  50  degrees  is  a 
strong  indication  of  a  bent  transition  state.38 

When  the  extrapolation  of  KDIE  is  carried  out  to  extremely  high 
temperatures  the  exponential  portion  of  equation  14  reduces  to  unity.  This 
leaves  only  the  preexponential  AH/AD  to  contribute  to  the  KDIE.  What  this 
preexponential  actually  represents  is  the  ratio  of  partition  functions  (eq.  17) 
which  at  high  temperature  reduce  to  v^L/v5L,  the  ratio  of  imaginary  frequencies 
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describing  the  motion  along  the  reaction  coordinate.22-32  At  infinite  temperature 
the  value  of  this  ratio  of  imaginary  frequencies  and  thus  the  value  of  the  KDIE 
should  in  theory  fall  between  unity  and  (p&/p&.)1/2  (the  ratio  of  the  reduced 

masses  along  the  reaction  path).32 

For  primary  KDIEs,  where  the  isotopic  atom  is  in  motion  along  the 
reaction  path,  the  infinite  mass  approximation  (see  equations  6,  7,  and  8)  can 
be  used  to  give 


where  m  is  the  mass  of  the  isotopic  atom.  Thus,  the  high  temperature  limiting 
value  for  a  primary  KDIE  should  be  (2)1/2  or  1.41.  However,  due  to 
anharmonicities  and  inaccuracies  in  the  infinite  mass  approximation  this  value  is 
closer  to  1.35.44 

In  secondary  KDIEs,  the  isotopic  atom  is  not  involved  in  motion  along  the 
reaction  path.  As  such  the  infinite  mass  approximation  cannot  be  used  to 
simplify  the  ratio  of  reduced  masses  (h$i/^l)1/2-  Because  of  the  small  size  of 
the  isotopic  atom  and  its  rather  minor  contribution  to  the  reduced  mass,  pJLis 
only  slightly  larger  than  p^L,  giving  a  value  of  (M&/p&.)1/2  very  close  to  unity. 
Thus,  at  infinite  temperature,  secondary  KDIEs  approach  unity. 

The  above  discussion  of  temperature  dependence  Is  an  oversimplification 
of  a  very  complex  subject.  Temperature  dependence  is  not  as  straightforward  as 
the  above  discussion  would  indicate.  However,  in  general,  one  can  conclude 
that  the  larger  the  value  of  [AEj£  for  a  particular  reaction  the  more  the  KDIE  for 

that  reaction  will  change  with  changes  in  temperature. 

CONCLUSIONS 

In  the  discussions  thus  far,  I  have  tried  to  focus  my  attentions  on  the  KDIE 
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arising  from  a  singla  isotopic  atom.  It  is  frequently  the  case,  however,  that 
reactions  are  performed  on  molecules  containing  more  that  one  isotopic  atom. 
This  gives  rise  to  the  possibility  of  compounded  isotope  effects.  Essentially, 
when  there  is  more  than  one  isotopic  atom,  the  effects  of  each  atom  multiply. 23 
This  of  course  assumes  that  there  is  only  one  isotopic  atom  contributing  a 
primary  effect  and  the  remaining  atoms  contribute  secondary  effects  or  that  all 
isotopic  atoms  contribute  secondary  effects. 

In  looking  back  at  the  calculation  of  the  secondary  KDIE  performed  by 
Streitwieser  et  al.,44  suppose  instead  of  using  one  isotopic  atom  they  had  used 
three.  The  resulting  calculations  would  involve  three  frequencies  from  each 
isotopic  atom.  The  resulting  KDIE  would  be  the  product  of  the  contributions  of  all 
nine  frequencies,  or  in  other  words  the  product  of  the  isotope  effects  for  each  of 
the  three  atoms.  If  the  isotope  effect  of  each  atom  were  the  same  as  the 
secondary  isotope  effect  for  the  single  atom  case  (i.e.,  1 .48),  then  the  resulting 
KDIE  would  be  (1 .48)3  or  3.24. 

The  compounding  of  isotope  effects  is  not  limited  to  secondary  effects.  It 
is  quite  possible  to  have  a  compounded  isotope  effect  where  one  atom 
contributes  a  primary  effect  and  one  or  more  atoms  contribute  a  secondary 
effect. 

The  presence  of  multiple  isotopic  substitution  also  complicates  the 
temperature  dependence  of  observed  KDIE.  Because  the  effects  of  the  various 
isotopic  atoms  may  exhibit  different  degrees  of  temperature  dependence,  the  net 
temperature  dependence  of  the  observed  KDIE  may  be  quite  complex. 

The  determination  of  whether  an  isotope  effect  is  primary  or  secondary  is 
relatively  straightforward  for  reactions  involving  a  single  isotopic  atom.  If  the 
observed  effect  is  greater  than  the  maximum  possible  secondary  effect,  then  the 
KDIE  can  unambiguously  be  assigned  to  a  primary  effect.  If  however,  the  KDIE 
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is  less  than  the  maximum  secondary  effect,  then  other  information  must  be  used 
to  determine  whether  the  effect  is  primary  or  secondary. 

The  evaluation  of  KOIE  becomes  an  order  of  magnitude  more  complex 
when  multiple  isotopic  substitution  is  involved.  Because  of  the  possibility  or 
large  compounded  secondary  effects,  only  the  largest  KDIE  can  be  confidently 
assigned  as  primary.  For  smaller  isotope  effects  it  may  be  possible  to  make 
clear  assignments  only  if  information  from  other  types  of  experiments  is  used. 


CRITICAL  EVALUATION  OF  ENERGETIC  MATERIAL  KDIE  STUDIES 
OVERVIEW 

Because  of  the  complexity  of  KDIE  theory,  approximations  have  been 
made  in  an  effort  to  simplify  its  application  to  experimental  data.  Unfortunately, 
KDIE  theory  has  frequently  been  oversimplified  to  the  point  where  it  becomes 
valid  only  under  idealized  circumstances.  The  resulting  application  of  this 
simplified  model  to  experimental  data  has  led  to  frequent  errors  in  the 
interpretation  of  KDIE.  26 

The  KDIE  studies  of  energetic  materials  epitomize  the  problems  in 
interpretation  which  have  plagued  isotope  effect  studies  since  their  inception. 
Publications  in  this  area  have  utilized  1.4146  (1.35  in  more  recent  papers)10  as  a 
dividing  line  between  primary  and  secondary  isotope  effects.  In  these  studies, 
experimental  KDIE  values  were  generally  designated  as  primary  if  they  were 
greater  than  1.41  (1.35)  and  as  secondary  if  they  were  less  than  1.41  (1.35). 
This  number  is  essentially  empirical  in  nature  and  possesses  only  limited 
theoretical  significance  with  respect  to  w*  ether  *n  effect  is  primary  or 
secondary.  Its  only  importance  in  KDIE  theory  stems  from  the  fact  that  at  infinite 
temperature  the  maximum  possible  primary  effect  is  approximately  1 J3 S.32-44  It 
is  quite  possible  to  have  compounded  secondary  effects  greater  than 
1.41(1.35),  and  it  is  also  altogether  possible  to  have  primary  effects  less  than 
1.41(1.35). 

Fortuitously,  the  1.41(1.35)  dividing  line  lies  within  the  range  of  what  can 
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be  expected  for  a  primary  effect  in  the  thermal  processes  of  energetic  materials. 
As  we  shall  see  in  the  following  discussion,  when  more  rigorous  guidelines  are 
applied  to  the  energetic  material  KDIE  data  the  assignments  made  using  the 
1 .41  (1 .35)  criterion  appear  in  most  cases  to  be  correct. 

The  application  of  KDIE  to  the  study  of  energetic  materials  has  been  at 
best  a  difficult  proposition.  In  addition  to  the  problems  associated  with  the 
oversimplification  of  KDIE  theory,  there  are  other  obstacles  to  obtaining  and 
interpreting  isotope  effect  data.  The  thermal  processes  of  energetic  materials 
(decomposition,  pyrolysis,  combustion,  and  detonation)  are  quite  complex. 
Their  mechanisms  consist  of  numerous  steps  and  may  involve  several  parallel 
pathways  each  with  its  own  rate  limiting  process.  The  precision  of  the 
techniques  used  to  study  these  thermal  processes  tends  to  be  poor,  which  in 
turn  leads  to  KDIE  values  with  relatively  large  uncertainties.  Often  kinetic  rate 
constants  cannot  be  measured  directly;  as  a  result,  ratios  of  other  parameters, 
which  are  hopefully  related  to  the  rate  constant,  must  be  used  to  obtain  isotope 
effect  values.  Finally,  the  very  nature  of  these  processes  necessitates  that 
experiments  be  performed  at  relatively  high  temperatures  (473  to  673K).  As  a 
consequence  of  this  and  other  factors,  KDIE  values  obtained  for  the  thermal 
processes  of  energetic  materials  are  smalt. 

TRINITROTOLUENE 

Shackelford  et  al.45  were  the  first  to  apply  the  technique  of  KDIE  to  the 
study  of  energetic  materials.  They  attempted  to  elucidate  the  rate  limiting 
process  in  the  thermal  decomposition  of  liquid  2,4,6*trinitrotolu&ne  (TNT).  Like 
most  energetic  organic  compounds,  TNT  exhibits  autocatalytic  behavior  upon 
decomposition.46-47  However,  unlike  the  nitramines  to  be  discussed  later,  TNT 
produces  relatively  few  gaseous  products  during  decomposition.48  The 
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decomposition  of  TNT  proceeds  initially  through  a  long  induction  period 
characterized  by  minimal  heat  evolution.  During  the  induction  phase  a  catalytic 
species  is  generated  which  at  a  certain  threshold  concentration  initiates  an 
exothermic  accelerator  phase.  The  accelerator  phase  is  followed,  after  the 
point  of  maximum  heat  release,  by  a  decay  phase  which  continues  until 
decomposition  is  complete.45*46 

Using  isothermal  differential  scanning  calorimetr  (IDSC),  Shackelford  et 
al.45  obtained  data  for  the  decomposition  of  TNT  and  its  a.a.ct'trideuteriomethyl 
analogue  (TNT-d3)  at  temperatures  ranging  from  518  to  542  K.  Applying  first 

order  kinetics  to  the  data  for  the  decay  phase,  they  obtained  rate  constants  for 
both  TNT  and  TNT-dj.  The  average  value  for  k^o  for  the  temperature  range 

studied  was  given  as  1 .35  ±  0.02.  However,  using  the  rate  constant  data  they 
supplied  for  each  temperature,  I  obtained  an  average  KDIE  for  the  decay  phase 
of  1 .38  ±  0.2C  at  the  95%  confidence  level.  No  explanation  can  be  offered  for 
the  disparity  in  the  manipulation  of  their  data.  (When  the  decay  phase  KDIE  is 
discussed  in  the  text  below,  both  isotope  effect  values  will  be  given.) 

Rate  constants  could  not  be  obtained  for  the  induction  phase  of  TNT 
decomposition.  However,  it  was  suggested  that  the  induction  time,  which  is  the 
time  between  the  beginning  of  the  DSC  experiment  and  the  onset  of  the 
accelerator  phase,  should  be  inversely  proportional  to  the  rate  constant  such 
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where  t)D  is  the  induction  time  for  TNT-d3  and  t,H  is  the  induction  time  for  TNT. 
Using  equation  19,  an  average  isotope  effect  for  the  induction  phase  of  1.66  ± 
0.20  was  obtained.45 

Shackelford  et  al.,  using  1.41  as  the  dividing  line  between  primary  and 
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secondary  KDIE,  stated  that  the  induction  phase  KP!E  of  1.66  was  clearly  a 
primary  isotope  effect.  Although  less  than  1 .41 .  they  attributed  the  decay  phase 
isotope  effect  of  1 .35  (1 .38)  to  a  diluted  primary  isotope  effect. 

The  following  chemical  evidence  was  used  to  support  the  above 
assignments.  Data  obtained  by  Dacons  et  al.40  and  Rogers49  indicated  that  the 
initial  decomposition  of  TNT  proceeds  via  the  oxidation  of  the  methyl  moiety. 
Depending  on  reaction  conditions  the  oxidation  of  TNT  may  proceed  from 

CH,  CH,OH  CHO  COOH 

O.N.A^NO,  no,  O.N  A  NO,  O^^^NO, 

—  w  —  —  oVW  <«» 

NO,  NO,  NO,  NO,  NO, 

TNT  1234 

2,4,6-trinitrobenzyl  alcohol  (1)  through  2,4,6-trimtrobenzyaldehyde  (2)  to  form 
2,4,6'trinitiobenzoic  acid  (3)  which  then  readily  decarboxylates  to  give  2,4,6* 
trinitrobenzene  (4)  (eq.  23),  or  equation  23  may  be  followed  through  the 
production  of  2  which  then  reacts  with  other  species  in  the  liquid  TNT  to 
produce  a  myriad  of  decomposition  products  called  "explosive  coke."40-49  The 
formation  of  1,2,  and  3  requires  the  rupture  of  C*H  bonds.  Although  the 
decarboxylation  of  3  necessitates  a  C-C  bond  rupture,  the  absence  of  a- 
hydrogens  precludes  any  significant  KDIE. 

When  small  amounts  of  hydroquinone,  which  possesses  a  labile 
hydroxyl  hydrogen,45  were  introduced  into  liquid  TNT  a  significant  decrease  in 
the  induction  time  was  observed.  These  results  appear  to  require  some  type  of 
bimolecular  hydrogen  transfer  between  hydroquinone  and  TNT,  and  thus  imply 
that  the  rate  controlling  process  during  the  induction  phase  (and  presumably 
during  the  entire  decomposition)  involves  a  hydrogen  transfer  between  two  TNT 
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molecules  (l.e.t  a  C*H  bond  rupture).  However,  one  cannot  entirely  rule  out  the 
possibility  that  the  hydroquinone  accelerates  a  radical  process  not  involving  a 
C-H  bond  rupture. 


X,  r  „ 

(Observed  KDIE-2°‘2°- 2°) 


NO, 


FIGURE  5 .  Possible  Isotope  Effects  in  TNTt^. 


For  a  TNT  molecule  with  deuterium  substitution  at  the  methyl  moiety,  an 
isotope  effect  should  be  observed  if  either  a  methyl  carbon-hydrogen  bonj  or 
the  benzylic  carbon-carbon  bond  were  broken  in  the  rate  determining  step.  In 
the  case  of  a  C-H  bond  rupture  there  would  be  two  secondary  effects,  due  to  the 
two  remaining  deuteriums  on  the  methyl  carbon,  in  addition  to  the  obvious 
primary  effect  (Fig.  5).  The  observed  isotope  effect  would  be  the  product  of  a 
primary  effect  compounded  by  two  secondary  effects  (i.e.,  1°*  2°*  2°).  As 
shown  in  Figure  5,  the  observed  isotope  effect  for  a  C-C  bond  rupture  would 
simply  be  the  product  of  the  secondary  effects  resulting  from  the  three  methyl 
deuteriums  (i.e.,  2°  •  2°  •  2°). 

The  TNT  decomposition  KDIE  of  1.35  (1.38)  arid  1.66  are  very  small  for 
compounded  primary  effects  when  considering  that  they  are  the  product  of  the 
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isotope  effects  from  three  deuterium  atoms.  Room  temperature  compounded 
secondary  effects  have  been  known  to  have  values  larger  than  the  KDIE 
observed  for  TNT.23*27*40*43  Consequently,  even  with  chemical  evidence 
supplied  by  Shackelford  et  al.,  the  absolute  assignment  of  the  TNT  KDIE  as 
primary  is  somewhat  tenuous. 

Fortunately,  the  evaluation  of  the  KDIE  in  TNT  decomposition  has  been 
significantly  clarified  by  additional  research  in  this  area.  Bulusu  and 
coworkers18*50  qualitatively  confirmed  the  work  of  Shackelford  et  al.46  by 
observing  isotope  effects  in  both  the  decomposition  and  detonation  of  TNT. 
More  importantly,  recent  electron  paramagnetic  resonance  (EPR)  studies51 
indicate  that  the  transition  state  for  TNT  decomposition  involves  a  bimolecular 
hydrogen  atom  transfer  between  the  methyl  moiety  of  one  TNT  molecule  and 
the  para  nitro-  group  of  a  second  TNT  molecule  (Fig.  6).  The  EPR  results,  which 
are  in  agreement  with  the  previous  hydroquinone  experiments,  provide 
reasonably  conclusive  proof  that  a  C-H  bond  rupture  is  indeed  the  rate 
controlling  process  in  TNT  decomposition. 


FIGURE  6.  Possible  Bimolecular  TNT  Transition  State. 


Interestingly,  the  ratio  of  preexponentials  (AH/AD)  for  the  decay  phase  of 
TNT  decomposition  is  approximately  one.  Assuming  the  decay  phase  KDIE  is 
temperature  dependent,  then  following  Kwart's36  criteria,  the  decomposition 
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should  proceed  through  a  linear  unsymmetric;  i  transition  state,  as  is  indicated 
by  EPR.  Unfortunately,  little  weight  can  be  givei.  to  this  observation  because  of 
the  large  uncertainties  and  undefined  temperature  dependence  of  the  TNT 
isotope  effect. 

Since  we  can  say  with  some  confidence  that  the  1.35  (1.38)  and  1.66 
values  are  compounded  primary  effects,  it  is  possible  to  make  estimates  of  the 
magnitude  of  the  primary  and  secondary  contributions  to  these  effects.  The 
thermolysis  of  4-methylene- 1 -pyrozolines  shows  an  average  secondary  effect  of 
1 .07  to  1 .08  per  a-deuterium  at  443  K.52  Assuming  this  reaction  to  have  a 
similar  transition  state  to  that  of  TNT  decomposition,  a  reasonable  estimate  for 
the  secondary  effect  would  he  1 .07  per  deuterium. 

In  the  few  studies  for  which  both  primary  and  secondary  effects  were 
independently  measured  for  the  same  reaction,  the  primary  isotope  effect  was 
always  significantly  larger  than  the  secondary  icotope  effect.23  62  For  example, 
Reference  23  gives  independently  measured  primary  and  secondary  Isotope 
effects  for  an  elimination  reaction  at  323  K  of  2.98  and  1.30,  respectively. 
Assuming  a  similar  relationship  between  the  primary  and  secondary  effects  for 
TNT  decomposition  and  that  the  secondary  effect  for  TNT  is  approximately  1 .07, 
a  corresponding  primary  effect  can  be  calculated  by  a  temperature 
extrapolation.  This  is  accomplished  by  first  using  equation  16  and  treating  the 
preexponential  ratios  as  one;  the  KDiE  of  1.30  at  298  K  is  extrapolated  to  the 
temperature  at  which  it  decreases  to  1 .07.  The  KDIE  of  2.98  at  298  K  is  then 
extrapolated  to  this  same  temperature,  giving  a  value  of  1.32.  Thus,  the 
estimated  primary  and  secondary  isotope  effects  for  TNT  decomposition  (518  to 
542  K)  per  deuterium  are  1.32  and  1.07,  respectively.  Using  these  estimates, 
one  would  obtain  a  value  for  the  compounded  primary  effect  in  TNT 
decomposition  of  1.51  (i.e.  1.32  •  1.07  •  1.07),  which  is  in  reasonable  agreement 
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with  both  the  observed  KDIE  values. 

It  must  be  pointed  out  that  by  simply  assuming  the  TNT  secondary  effect 
to  be  1 .07  and  then  determining  what  magnitude  of  primary  effect  would  give 
the  observed  KDIE  values,  a  value  of  approximately  1 .32  for  the  primary  effect 
would  be  obtained. 

The  magnitude  of  the  contribution  of  a  deuterium  atom  to  an  observed 
KDIE  will  certainly  depend  on  the  configuration  of  the  transition  state  and  clearly 
reactions  involving  different  compounds  will  have  different  transition  states. 
However,  there  are  some  similarities  between  the  thermal  reactions  of 
nitroaromatics  and  nitramines.  Although  crude,  the  TNT  KDIE  values  of  1 .32  for 
a  primary  contribution  and  1 .07  for  a  secondary  contribution  provide  the  only 
reasonable  estimate  for  isotope  effect  values  in  the  thermal  processes  of 
energetic  materials  at  elevated  temperatures.  Consequently,  these  values  will 
be  used  to  help  evaluate  KDIEs  in  the  following  discussions. 

TRIAMINOTRINITROBENZENE 

Rogers  et  al,53  extended  the  isotope  effect  studies  of  energetic  materials 
to  the  solid  state  decomposition  of  1 , 3, 5-triamino-2, 4, 6-trinitrobenzene  (TATB). 
Performing  IDSC  experiments  on  TATB  and  TATB -cf9  similar  to  those  done  for 
TNT  over  a  temperature  range  of  570-655  K,  they  obtained  a  value  for  kH/k0  of 
approximately  1.6.  Using  the  1.41  citerion,45  they  attributed  this  to  a  primary 
effect. 

Because  of  the  lack  of  information  given  in  their  paper  it  is  difficult  to 
assess  the  validity  of  the  conclusions  made  by  Rogers  et  al.53  The  isotope  effect 
in  TATB  is  due  to  changes  in  the  force  constants  of  N-H  bonds,  as  opposed  to 
C-H  bonds  in  the  case  of  TNT.  N-H  stretching  frequencies  are  slightly  larger 
than  the  corresponding  C-H  frequencies,54  so  one  would  expect  the  isotope 
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effects  for  N-H  bonds  in  general  to  be  somewhat  larger  than  those  for  C-H 
bonds.  However,  due  to  the  fact  that  the  TATB  decomposition  occurs 
approximately  100  degrees  higher  the  TNT  decomposition,  our  reference 
values  of  1 .32  for  primary  and  1 .07  for  secondary  should  still  be  approximately 
vaHd.  This  of  course  assumes  that  the  TATB  KDIE  decreases  with  increasing 
temperature.  Thus,  if  the  numerous  approximations  are  reasonable,  then  a 
primary  KDIE  in  TATB  decomposition  stiouid  be  around  1.41  (i.e.,  1°  •  2°  •  1.32 
•  1.07).  Considering  the  precision  of  IDSC  data  and  the  approximations  made, 
this  is  in  relatively  good  agreement  with  experiment. 

The  results  of  Rogers  et  al.54  are  important  in  that  they  were  the  first  to 
demonstrate  a  significant  isotope  effect  in  the  solid  state.  The  solid  state  TATB 
experiments  differ  from  the  liquid  TNT  experiments  In  that  the  TATB  is  much 
more  restricted  in  its  motions  and  thus  more  limited  ir,  the  possible  transition 
state  configurations  it  can  assume,  due  to  the  strong  intermolecular  forces  in  the 
solid  state.  Therefore,  it  is  somewhat  surprising  that  a  KDIE  of  this  magnitude 
was  observed  for  TATB  decomposition. 

HMX  AND  RDX 

Recently,  KDIE  have  been  usao  to  study  the  rate  limiting  processes  in  the 
decomposition,17'20  deflagration,19  and  combustion16-17-56  of  HMX  and  RDX. 
The  evaluation  of  the  isotope  effects  and  assignment  of  the  rate  limiting 
processes  for  these  thermal  events  was  made  using  almost  exclusively  the 
empirical  criteria  of  Shackelford  et  al.  (i.e.  1.35). 19  Unlike  the  TNT  case, 
however,  there  is  lack  of  clear  experimental  evidence  to  corroborate  the 
determinations  made  from  the  isotope  effect  data.  Consequently,  conclusions 
made  from  KDIEs,  which  will  be  discussed  below,  must  be  treated  with  care 
until  other  methods  can  be  found  to  verify  these  results. 
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HMX  and  RDX  are  cyclic  nitramines  which  differ  by  only  one  methylene 
nitramine  unit.  Consequently,  the  possible  isotope  effects  for  HMX  and  ROX  are 
almost  identical.  For  deuterium  substitution  at  the  methylene  carbons,  isotope 
effects  should  be  observed  if  either  a  C-H,  C-N,  or  N-N  bond  were  broken  in  the 
rate  limiting  process.  A  C-H  bond  rupture  would  result  in  a  compounded 
primary  effect  due  to  a  primary  contribution  from  the  deuterium  in  motion  along 
the  reaction  path  and  a  secondary  contribution  from  the  deuterium  alpha  to  the 
reaction  path  (Fig.  7).  Using  the  TNT-based  values,  the  estimated  magnitude 
for  a  KDIE  from  a  C-H  bond  rupture  would  be  1.32  •  1.07  ■  1.41  (by  sheer 
coincidence  this  haopens  to  be  the  same  as  the  empirical  value  originated  by 
Shackelford  et  al.46).  In  theory,  the  four  deuteriums  gamma  to  the  reaction  path 
could  contribute  to  the  compounded  primary  effect.  However,  considering  the 
small  size  of  a-secondary  effects  under  these  conditions,  It  should  be 
reasonable  to  neglect  the  secondary  effects  of  the  ^-deuteriums. 

A  C-N  bond  rupture  in  the  rate  limiting  process  would  result  in  a 
compounded  secondary  effect  made  up  of  the  contributions  from  the  two  a- 
deuteriums  and  the  two  ^  deuteriums  (Fig.  7).  Depending  on  the  orientation  of 
the  activated  complex,  the  isotopic  effects  of  (3-deuteriums  (secondary  isotope 
effects  of  the  second  kind)  can  be  as  large  as  the  effects  of  a-deuteriums 
(secondary  isotope  effects  of  the  first  kind).  Under  most  conditions,  however, 
the  effects  of  (3-deuteriums  are  smaller.  Assuming  this  to  be  the  case,  an 
acceptable  estimate  for  the  contribution  of  the  (3-deuteriums  to  the  compounded 
secondary  KDIE  would  be  1 .02  per  deuterium.23-40  43  Using  this  approximation 
and  the  TNT  based  values,  the  estimated  magnitude  for  the  effect  for  a  C-N 
bond  rupture  would  be  1 .07  •  1 .07  •  1 .02  •  1 .02  -  1.19.  For  this  estimate  I  have 
chosen  to  neglect  the  possible  contributions  of  deuteriums  gamma  and  delta  to 
the  reaction  path. 
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There  is  sufficient  evidence  both  theoretical76  and  experimental15'10  20  to 
conclude  that  the  N-N  bond  is  the  weakest  covalent  bond  in  HMX  and  RDX. 
Furthermore,  these  studies  indicate  that  the  N-N  bond  rupture  is  relatively  facile 
at  the  temperatures  of  nitramine  decomposition.56  Consequently,  it  was 
concluded  that  the  N-N  bond  rupture  was  unlikely  to  be  the  rate  limiting 
process.  For  the  sake  of  this  research  we  will  assume  this  conclusion  to  be 
correct. 
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FIGURE  7.  Possible  Isotope  Effects  in  RDX-de  and  HMX-cfe. 


KDIE  in  HMX  Decomposition 

HMX  decomposition,  as  measured  by  IDSC,  possesses  three  phases 
similar  to  what  is  seen  in  TNT  decomposition.  Shackelford  et  al.19  attributed  the 
induction  phase,  the  exothermic  acceleratory  phase,  and  the  decay  phase  to 
physical  states  of  HMX:  solid,  nixed  melt,  and  liquid,  respectively.  The  KDIEs 
for  each  of  these  phases  are  listed  in  Table  I. 

There  are  two  KDIE  values  listed  for  each  phase  of  HMX  decomposition. 
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The  values  given  for  HMX  and  HMX-da  prepared  under  different  conditions 
have  to  be  questioned  somewhat  because  it  is  not  Known  how  much  of  the  rate 
differences  are  due  to  differences  in  sample  history.  Therefore,  the  KDIE  values 
for  the  identically  prepared  HMX  and  HMX-dg  should  be  given  more  credence 

even  though  they  represent  a  narrower  temperature  range  and  smaller  number 
of  replicate  measurements. 


TABLE  I 


Deuterium  Isotope  Effects  in  HMX  Decomposition 


Method 

Temp.  Range  (K) 

Physical  State 

Ratio 

KDIE 

Ref. 

IDSC 

551-553 

solid 

2.21  ±0.18* 

19 

IDSC 

551 

solid 

t»c/ln 

1.74  ±  0.22b 

19 

IDSC 

551-553 

mixed  melt 

W<D 

0.85  ±  0.22* 

19 

IDSC 

551 

mixed  melt 

W<D 

1.05  ±  0.18b 

19 

IDSC 

551-553 

liquid 

»VKd 

1.13  ±0.08* 

19 

IDSC 

551 

liquid 

W<D 

1.28  ±  0.2lb 

19 

TGA 

508-555 

solid  and  melt 

W<D 

2.07±0.11* 

20 

TGA 

508-555 

solid  and  melt 

1.91±  0.17* 

20 

•  HMX  and  HMX-d8  were  prepared  by  different  methods.  b  HMX  and 
HMX-d0  were  prepared  by  identical  methods. 


Shackelford  et  al.,19  using  their  empirical  criteria,  attributed  the  induction 
phase  KDIEs,  2.21  ±  0.18  and  1.74  ±  0.22,  to  a  rate  controlling  C-H  bond 
rupture.  To  determine  whether  these  values  are  indeed  due  to  a  compounded 
primary  effect,  one  needs  to  look  at  what  it  would  take  for  these  values  to  be  the 
result  of  a  compounded  secondary  effect.  Compounded  secondary  effects  have 
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bean  Known  to  be  larger  than  these  values  but  not  at  the  temperatures  of  HMX 
decomposition.  Even  if  one  assumed  that  the  secondary  effects  for  the  (5- 
deuteriums  were  as  large  as  those  for  the  a-deuteriums  (i.e.,  1.07),  one  would 
still  obtain  only  a  compounded  secondary  effect  of  (1.07)4  or  1.31,  which  is  far 
below  the  observed  value.  To  arrive  at  a  compounded  secondary  effect  in  the 
range  of  1.74  would  require,  for  example,  effects  of  1.20  for  each  a-deuterium 
and  1.10  for  each  {J-deuterlum.  While  not  impossible,  secondary  effects  of 
these  magnitudes  are  unlikely  at  the  elevated  temperatures  of  HMX 
decomposition.  Therefore  by  default,  it  would  seem  that  the  assessment  of  a 
rate  controlling  C-H  bond  rupture  is  reasonable. 

The  values  given  for  the  KDIE  in  the  mixed  melt  were  0.85  ±  0.22  for  the 
HMX  and  HMX-cf8  prepared  by  different  methods  and  1.05  ±  0.18  for  the 

identically  synthesized.  Shackelford  et  al.,18  placing  more  significance  on  the 
0.85  value,  stated  that  the  mixed  melt  exhibited  an  inverse  isotope  effect  which 
they  attributed  to  differences  in  C-H  and  C-D  bond  contractions  during 
liquefaction.  Because  their  explanation  places  a  great  deal  of  significance  on 
the  physical  state  of  the  HMX,  it  would  seem  that  the  value  for  the  identically 
synthesized  HMX  and  HMX-da  would  be  a  more  valid  measure  of  this  type  of 

isotope  effect.  The  rate  differences  giving  this  "inverse”  isotope  effoct  of  0.85 
are  likely  the  result  of  differences  in  impurities,  crystal  size,  defects,  etc. 
between  the  HMX  and  HMX-d8  and  not  the  result  of  deuterium  substitution. 

Given  the  conflicting  data  and  the  large  uncertainties  in  both  KOIEs  for  the 
mixed  melt,  the  no  conclusive  evidence  of  a  KDIE  in  the  mixed  melt,  inverse  or 
normal. 

The  KDIE  values  for  the  liquid  HMX  (1.13  ±  0.08  and  1.28  ±  0.21)  were 
attributed  by  Shackelford  et  al.  to  a  rate  controlling  C-N  bond  rupture.  However, 
isotope  effects  of  these  magnitudes  could  quite  possibly  be  either  compounded 
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primary  or  compounded  secondary  effects.  Without  other  information  it  is 
impossible  to  confirm  their  assignment. 

Bulusu  et  al.20  employed  thermogravimetnc  analysis  (TGA)  to  study 
isotope  effects  in  HMX  decomposition.  Because  nitramine  decomposition 
proceeds  to  gaseous  products  early  on,  thermogravimetric  analysis,  which 
measures  sample  weight  loss  with  time,  is  a  viable  method  for  measuring  rate 
constants  and  thus  isotope  effects.  The  KDIE  values  obtained  for  HMX 
decomposition  using  TGA  are  listed  in  Table  I.  Interestingly,  Bulusu  et  al.  did 
not  observe  the  complex  decomposition  behavior  that  was  seen  in  IDSC 
experiments.  They  obtained  essentially  one  KDIE  for  the  entire  decomposition 
experiment  (2.07  ±  0.11  determined  from  the  ratio  of  rates  and  1.S  ±  0.17 
determined  from  the  ratio  of  induction  times).  Using  the  empirical  1.35 
criterion,19  Bulusu  et  al.  attributed  the  isotope  effect  to  a  rate  limiting  C-H  bond 
rupture.  Following  the  logic  used  to  explain  the  KDIE  in  the  IDSC  induction 
phase,  the  isotope  effects  measured  by  TGA  by  default  appear  to  be  a  primary 
effect  compounded  by  secondary  effects. 

One  interesting  observation  made  by  Bulusu  et  al 20  was  the  apparent 
temperature  independence  of  the  HMX  KDIE  over  a  temperature  range  of  508- 
555  K.  From  this  observation  they  attempted  to  invoke  Kwart's38  criteria  to 
suggest  that  the  HMX  decomposition  went  by  a  non-linear  transition  state.  This 
observation,  however,  is  little  more  than  conjecture.  Data  an  order  of 
magnitude  more  precise  than  what  is  obtained  from  either  the  TGA  or  IDSC 
experiments  would  be  necessary  before  any  justifiable  conclusions  about  the 
configuration  of  the  transition  state  in  HMX  decomposition  could  be  made. 

KDIE  in  RDX  Decomposition 

The  KDIE  values  obtained  for  RDX  decomposition  -ire  given  in  Table  II. 
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The  isotope  effect  values  for  IDSC  and  TQA  are  in  the  same  range  even  though 
the  TQA  measured  the  KDIE  for  the  solid  state  decomposition  while  IDSC 
measured  KDIE  for  decompos’tion  in  the  liquid.  Using  the  1.35  criterion,19  ail 
three  KDIE  for  RDX  decomposition  were  designated  as  primary.17*20  Although 
the  TQA  ualue  of  1.50  is  smaller  than  the  IDSC  values,  it  would  still  require  a 
relatively  large  secondary  effect  from  all  four  deuteriums  for  it  to  be  the  result  of 
a  compounded  secondary  effect.  Therefore,  as  was  the  case  of  the  HMX 
induction  phase  because  it  seems  unlikely  that  the  effects  can  be  due  to  a 
compounded  secondary  effect,  the  KDIE  for  RDX  decomposition  appears  to  be 
the  result  of  a  rate  controlling  C-H  bond  rupture. 


TABLE  II 

Deuterium  Isotope  Effects  in  RDX  Decomposition 


Method 

Temp.  Range  (K) 

Physical  State 

Ratio 

KDIE 

Ref. 

IDSC 

505 

liquid 

M<d 

2.06±  1.07*b 

17 

IDSC 

510 

liquid 

Kh^o 

1.74  ±  0.1 7* 

17 

TQA 

472-489 

solid 

W<d 

1 .50  ±  0.06c 

20 

*  This  is  a  preliminary  value  based  on  six  measurements  each  for 
RDX  and  RDX-dfl.  b  The  RDX  and  RDX-de  were  prepared  by  identical 
methods.  c  The  RDX  and  RDX-d9  were  prepared  by  different  methods. 


KDIE  in  HMX  and  RDX  Combustion 

The  burn  rates  of  pressed  pellets  uf  HMX  and  HMX-ds  were  measured  at 
three  different  pressures.  The  ratio  of  the  burn  rates  (r^r^,  which  was  assumed 
to  be  equivalent  to  the  ratio  of  rate  constants,  at  each  pressure  are  given  in 
Table  III.  Shackelford  et  at.,16  using  the  1.35  dividing  line,  designated  the  1.37 
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±  0.24  and  1.61  ±  0.11  at  500  and  1000  psig  respectively  as  primary  effects. 
The  value  at  1 500  psig  of  1 .24  ±  0. 1 3  was  assigned  to  a  secondary  effect. 

HMX  combustion  normally  proceeds  through  an  exothermic 
decomposition  producing  gaseous  products  which  further  react  in  the  flame. 
The  change  in  the  isotope  effect  with  pressure  was  attributed  to  a  change  in  the 
physical  state  where  a  majority  of  the  HMX  decomposition  occurred. 
Shackelford  et  al.16  speculated  that  as  pressure  increased,  the  HMX 
decomposition  went  from  the  solid  state  to  the  liquid.  While  their  explanation 
for  the  variation  in  KOIE  values  is  possible,  it  is  much  more  likely  that  the 
differences  between  pressures  are  simply  due  to  statistical  fluctuations.57  If  this 
is  correct,  then  a  more  reasonable  value  for  the  isotope  effects  in  HMX 
combustion  would  be  an  average  of  the  KDIE  for  all  three  pressures  (i.e.,  1.40  ± 
0.11). 


TABLE  III 

Deuterium  Isotope  Effects  in  HMX  and  RDX  Combustion 


Compound 

Pressure  (psig) 

KDIE  (rH/rD) 

Ref. 

HMX 

500 

1.37±0.24 

16 

HMX 

1000 

1.61  ±0.11 

16 

HMX 

1500 

1.2410.13 

16 

RDX 

500 

1.3710.16 

17 

RDX 

1000 

1.4610.23 

17 
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The  bum  rates  of  pressed  pellets  of  RDX  and  RDX-de  were  also 
measured  at  500  and  1000  psig.  The  ratio  of  these  burn  rates  at  each  pressure 
is  listed  in  Table  III.  Shackelford  et  al.,17  using  the  same  logic  applied  to  HMX, 
designated  these  effects  as  primary. 

The  isotope  effects  in  HMX  and  RDX  combustion  are  larger  than  the 
estimated  magnitude  for  the  compounded  effect  recuiting  from  a  C-N  bond 
rupture.  However,  given  the  uncertainties  in  both  the  KDIE  values  and  the 
compounded  secondary  estimate,  this  is  far  from  conclusive  proof  that  the 
combustion  isotope  effects  are  the  result  of  a  C-H  bond  rupture  in  the  rate 
limiting  step. 

If  one  assumes  that  the  rate  of  decomposition  of  HMX  and  RDX  in  the 
condensed  phase  controls  the  rate  of  combustion,  then  the  KDIE  results  for 
decomposition  can  be  used  to  help  interpret  the  isotope  effects  for  HMX  and 
RDX  combustion.  The  average  HMX  KDIE  of  1 .40  ±  0.1 1  and  the  average  RDX 
KDIE  of  1 .42  ±  0.1 5  are  smaller  than  the  compounded  primary  isotope  effects  in 
HMX  and  RDX  decomposition.  Much  of  this  difference  may  simply  be  due  to  the 
fact  that  the  temperature  of  the  burning  surface,  where  nitramine  decomposition 
takes  place,  is  approximately  100  to  200  degrees  higher  than  the  temperatures 
studied  for  HMX  and  RDX  decomposition.  Assuming  even  a  small  temperature 
dependence  of  the  KDIE  for  HMX  and  RDX  decomposition,  a  change  in 
temperature  of  this  magnitude  could  account  for  the  differences  in  isotope  effect 
values  between  decomposition  and  combustion.  If  this  is  truly  the  case,  then  it 
would  indicate  that  the  KDIEs  for  HMX  and  RDX  combustion  are  more  than 
likely  the  result  of  a  rate  controlling  C-H  bond  rupture. 

The  true  relationship  between  burn  rate  and  the  kinetic  rate  constant  of 
the  slowest  step  has  yet  to  be  clearly  addressed.  Shackelford  and 
coworkers16'17  simply  assumed  the  two  quantities  were  proportional  to  each 
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other.  While  reasonable,  this  assumption  needs  to  be  investigated  in  more 
detail. 

CONCLUSIONS 

The  poor  precision  of  the  KOIE  for  the  decomposition  and  combustion  of 
the  nitramines  (HMX  and  RDX)  and  the  lack  of  other  chemical  evidence  would 
make  many  chemists  shy  away  from  making  many  conclusions  concerning  the 
possible  rate  limiting  process.  One  must  understand  the  difficulties  associated 
with  these  types  of  experiments  and  accept  the  above  data  as  the  best  available 
at  the  present  time.  One  can  say  with  confidence  that  the  isotope  effects  are 
real.  However,  to  go  beyond  this  requires  that  the  assumptions  and 
approximations  made  in  the  above  discussion  be  valid.  If  they  are,  then  it  would 
appear  that  under  most  circumstances,  a  C-H  bond  rupture  is  the  rate 
controlling  process  in  HMX  and  RDX  decomposition  and  combustion. 
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EXPERIMENTAL 


SYNTHESIZED  COMPOUNDS 

This  experimental  section  describes  the  synthesis  of  HMX,  RDX,  and 
their  fully  deuteriated  analogs  HMX-d8  and  RDX-d8.  The  synthesis  of  HMX  and 

RDX  are  described  in  varying  degrees  of  detail  in  the  literature.1 61 9  58  62 
However,  for  the  sake  of  completeness  and  at  the  risk  of  repeating  information 
already  contained  in  the  literature,  I  will  endeavor  to  give  a  full  account  of  the 
synthesis  of  HMX,  HMX-df8,  RDX,  RDX-c/e. 

To  ensure  that  differences  in  experimental  results  were  not  due  to 
impurities  or  sample  history,  the  deuteriated  forms  of  HMX  and  RDX  were 
synthesized  in  the  same  manner  as  the  protiated  forms.  This  was  accomplished 
by  simply  substituting  paraformaldehyde-^?  for  paraformaldehyde  in  the 
synthesis  of  hexamethylenetetramine  to  give  hexamethylenetetramine-d12.  The 

deuteriated  hexamethylenetetramine  was  then  used  in  the  synthesis  of  either 
HMX-da  or  RDX-d6 

All  chemicals  were  reagent  grade  or  better  and  were  used  without 
purification.  The  paraformaldehyde-^?  was  obtained  from  MSD  Isotopes. 

Burdick  and  Jackson  distilled-in-glass  acetone  was  used  to  convert  HMX  fros.i 
the  a-  to  the  p-polymorph.  A  JEOL  FX-90Q  90MHz  FT-NMR  was  used  to  obtain 
all  ’H  spectra.  Melting  point  data  were  obtained  using  a  Thomas-Hoover 
Unimelt  capillary  melting  apparatus  and  unsealed  capillary  tubes.  All  melting 
points  are  uncorrected.  An  Olympus  model  AHMT  optical  microscope  was  used 
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to  obtain  optical  micrographs  of  the  a  and  p  forms  of  HMX. 

Synthesis  of  Anhydrous  HN03 

Because  98*100%  nitric  add  is  an  essential  reagent  in  the  final  nitration 
steps  of  both  HMX  and  RDX,  it  is  first  necessary  to  describe  its  preparation.63 
Concentrated  sulfuric  acid  (1000  mL)  was  added  to  a  two-liter  three-neck 
round-bottom  flask.  The  sulfuric  acid  was  cooled  in  a  dry  ice/acetone  bath  as 
600  g  (5.95  mol)  of  KN03  was  slowly  added  with  stirring  from  an  overhead 

mechanical  stirrer.  Upon  complete  dissolution  of  the  potassium  nitrate,  the 
mechanical  stirrer  was  removed,  and  a  still  head,  thermometer,  condenser,  bent 
adapter,  and  500-mL  round-bottom  flask  were  attached  in  preparation  for 
vacuum  distillation.  The  500-mL  collection  flask  was  placed  in  a  dry  ice/ 
acetone  bath,  and  a  methanol-ethylene  glycol  mixture  cooled  to  -30  °C  was 
circulated  through  the  condenser.  The  apparatus  was  then  slowly  evacuated  to 
a  pressure  of  46-47  torr.  The  reaction  flask  was  allowed  gradually  to  warm  to 
room  temperature  as  the  HN03  distilled  off.  After  HN03  bubble  evolution 

slowed,  the  solution  was  slowly  heated  to  60  °C.  The  distillation  was  stopped 
as  soon  as  any  evidence  of  red  N02  was  observed.  Normally  this  was  seen  as 
a  slight  yellow  tint  to  the  condensate  in  the  apparatus.  The  frozen  HN03  was 
allowed  to  liquefy  slowly  in  a  refrigerator.  The  clear  liquid  HNOa  was  then 
transferred  to  Teflon  bottles  under  a  nitrogen  atmosphere  and  stored  in  a 
refrigerator  to  prevent  decomposition. 

Yield  362  g  (97%). 

CAUTION:  HMX  and  RDX  are  both  high  explosives  and  as  such  great  care 
must  be  exercised  at  all  times  during  the  synthesis  of  these  compounds.  HMX 
is  1 .6  times  more  powerful  than  TNT  and  is  second  only  to  hexanitrobenzene  in 

47 


explosive  force.  Teflon  (not  Teflon  coated)  utensils  must  be  used  in  place  of 
metal  utensils  during  synthesis  to  preclude  accidental  spark  ignition.  In 
addition,  all  glassware  should  be  cleaned  thoroughly  and  be  free  of  scratches. 

tvm.  f  •-  HMX 

synthesis  of  HMX  and  HMX-d8  follows  a  four  step  procedure 
diagramed  in  Figure  8  and  outlined  in  detail  in  the  text  below. 
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FIGURE  8.  HMX  Synthesis. 

Hexamethylenetetramine  (HMTA1.16-58-60  To  11.58  g  (0.386  mol)  of 
paraformaldehyde  slurried  in  20  mL  of  distilled  water  was  added  20  mL  of 
concentrated  NH4OH  (28-30%  NH3)  over  a  20-30  minute  period  while 
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maintaining  the  reaction  temperature  at  25-35  °C.  After  the  addition  of  the 
ammonium  hydroxide  was  complete,  the  reaction  solution  was  heated  in  an  oil 
bath  at  35-40  °C  for  24  hours.  Water  was  removed  from  the  reaction  mixture 
under  vacuum  until  a  wet  solid  remained.  The  moist  solid  was  dissolved  in  200 
mL  of  methylene  chloride  giving  two  iayers.  The  CH2CI2  layer  was  removed 
and  the  remaining  aqueous  layer  was  extracted  with  4X50  mL  CH2Cl2.  The 
combined  methylene  chloride  extracts  were  dried  over  1  g  of  MgS04.  The 
aqueous  layer  was  then  added  to  100  mL  of  CH2CI2  and  dried  over  5  g  of 
MgS04.  The  methylene  chloride  solutions  were  filtered  and  the  MgS04  washed 
several  times  with  CH2CI2.  Following  filtration,  the  methylene  chloride  was 

removed  under  vacuum,  leaving  a  clean  white  powder. 

Yield  8.67  g  (96%),  m.p.  255-260  °C  (decomposes),  ’H  NMR  (CDCI3)  8 
4.66  (-CH2-). 

Hexamethylenetetramine-^  (HMTA-d12T60  The  HMTA-d12  was 
synthesized  using  15.00  g  (0.468  mol)  of  paraformaldehyde-d2  in  the  same 
manner  as  the  HMTA  described  above. 

Yield  1 1 .33  g  (95%),  m.p.  258-265  °C  (decomposes). 

1  ■5-dlacetvl-oclahvdro-3.7-endo-methvlene-1 .3.5.7-tetrazocine 
fDAPTI.80,81  To  8.65  g  (0.06  mol)  of  HMTA  was  added  3.51  g  (0.05  mol)  of 
NH4OAc  and  1.6  mL  of  distilled  water.  An  overhead  mechanical  stirrer  was 

used  to  mix  the  wet  solid  while  16.6  g  of  (0.016  mol)  of  acetic  anhydride  was 
added  over  a  one  hour  period.  During  the  addition  of  the  acetic  anhydride  the 
reaction  temperature  was  maintained  at  5-10  °C.  After  completion  of  the  acetic 
anhydride  addition,  the  stirred  solution  was  heated  at  35-40  °C  for  22  hours. 
The  reaction  yielded  approximately  25  mL  of  a  clear  viscous  liquid.  Half  of  this 


49 


liquid  was  used  for  the  synthesis  of  OAON  without  isolation  of  the  DAPT. 

1 .5-diacatvl-octahvdro-3.7-encfc>-methvl9ne-1 .3.5.7-tetrazoclna-cf10 
(DAPT-d10).  The  DAPT-d10  was  synthesized  using  9.39  g  (0.06  mol)  of  HMTA- 
d12.  The  remainder  of  the  synthesis  was  identical  to  the  synthesis  of  DAPT 
described  above. 

1  ■&-diacatvl-octflhvdro-3.7-dinltro«1 .3.5.7-tetrazocine_  fDAQNL,66°61  A 
total  of  16.98  g  (0.168  mot)  of  KN03  was  slowly  add9d  to  55  mL  of  concentrated 
H2S04  while  maintaining  25-30  °C.  Half  the  DAPT  solution  from  above  (-12.5 
mL)  was  then  added  to  the  acid  solution  at  a  rate  such  that  the  reaction 
temperature  remained  at  25-30  °C.  During  and  after  the  addition  of  the  DAPT 
solution,  vigorous  stirring  from  an  overhead  mechanical  stirrer  was  required  to 
effect  complete  mixing  of  the  two-phase  system.  When  the  addition  of  DAPT 
was  complete,  the  solution  was  stirred  for  an  additional  hour  at  25-35  °C.  The 
reaction  was  quenched  into  500  g  of  ice  water  and  then  placed  into  a 
refrigerator  where  the  DADN  product  crystallized  overnight.  The  product  was 
suction  filtered,  washed  with  water,  and  air  dried. 

Overall  yield  for  DAPT  and  DADN  7.44  g  (83%),  m.p.  264-265  °C.  ’H 
NMR  (DMSO-dy  of  DADN  gave  5  2.22  (-CH3),  8  5.51  (-CH2-). 

1  .S»diacetyl-octahydro-3.7-dinitro-1 .3.5.7-tatrazocine-d9  (DADN:.<fri.60 
The  DADNd8  was  synthesized  using  half  of  the  DAPT-  *10  solution  (-12.5  mL) 

in  the  same  manner  as  the  synthesis  of  DADN  described  above. 

Overall  yield  for  DAPT-d10  and  DADN-d8  7.46  g  (80%),  m.p.  264-265  °C, 

’H  NMR  confirmed  99%  isotopic  purity  calculated  from  the  acetyl  methyl 
protons. 

Qctahvdro-1 .3.5.7-tetranitro-l  .3.5.7-tetrazocine  fHMX).16-60'61  In  the 
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synthesis  of  HMX  the  configuration  of  the  reaction  apparatus  can  significantly 
affect  product  yield.  We  chose  to  use  a  relatively  large  reaction  flask  (500*mL 
three-neck  round-bottom)  to  allow  for  a  greater  volume  above  the  reaction 
mixture.  A  gas  bubbler  filled  with  high  vacuum  silicone  oil  was  used  to  limit  the 
venting  of  N02  and  N2Os  produced  during  the  reaction  and  prevent  the  inflow  of 

atmospheric  moisture.  The  conceptual  basis  of  the  apparatus  configuration  was 
to  maintain  a  positive  pressure  of  N02  and  N206,  the  actual  nitrating  agent,  in 

the  reaction  flask.  As  with  the  DADN  synthesis,  vigorous  stirring  was  employed 
by  means  of  an  overhead  mechanical  stirrer.  A  small  amount  of  anhydrous 
nitric  acid  was  used  to  lubricate  the  shaft  of  the  mechanical  stirrer. 

To  78  g  of  polyphosphoric  acid  was  added  35.3  g  of  anhydrous  HN03. 

These  acids  were  mixed  until  they  became  a  homogeneous  viscous  liquid.  A 
total  of  6.40  g  (0.022  mol)  of  DAON  was  added  rapidly  to  the  add  mixture.  The 
reaction  was  then  heated  at  60  °C  with  stirring.  As  the  reaction  progressed,  the 
atmosphere  above  the  solution  gradually  changed  from  a  light  red  initially  to  a 
dark  red-brown  by  the  end  of  the  reaction.  After  one  hour  the  reaction  was 
allowed  to  cool  to  40  °C  and  then  quenched  by  pouring  on  30C  g  of  finely 
crushed  Ice.  The  resulting  mixture  was  then  filtered,  washed  thoroughly  with 
water,  and  air  dried.  The  product  consisted  of  fine  white  needle-like  crystals 
(Fig.  9),  which  are  indicative  of  the  a  form  of  HMX.  The  (J-polymorph  of  HMX 
was  obtained  by  dissolving  the  a-HMX  in  acetone  and  then  gradually  removing 
the  solvent  with  vacuum.  The  conversion  to  p-HMX  was  confirmed  using  optical 
microscopy  (Fig.  10). 

Yield  p-HMX  4.98  g  (76%),  m.p.  279-280  °C,  ’H  NMR  (DMSO-de)  6  6.02 
(-CH2-).  The  DADN  methylene  peak  at  5  5.52  was  absent  from  ihe  NMR 
spectrum,  indicating  that  nitration  was  greater  than  98%  complete. 


51 


100  (UK 


Qctahvdro-1 .3.5.7-tetranitro-l .3.5.7-tetrazocine-d0  (HMX-d^.1680-61  The 
HMX*dt  was  synthesized  using  6.56  g  of  DADN-d6  in  an  identical  manner  to 
HMX  described  above.  Yield  p-HMX-d8  5.24  g  (e0%),  m.p.  279-280  °C,  ]H 
NMR  showed  small  peaks  at  5  5.52  and  6  6.02  corresponding  to  approximately 
0.5%  DADN-d^  and  the  equivalent  of  0.5%  HMX-h8  respectively. 

Synthesis  of  RDX 

The  synthesis  of  RDX  begins  with  the  preparation  of  HMTA  as  does  the 
synthesis  of  HMX.  However,  in  contrast  to  HMX,  the  preparation  of  RDX 
involves  the  conversion  of  HMTA  to  1,3,5-triacetyl-hexahydro-1,3,5-triazine 
(TRAT).  Furthermore,  the  final  nitration  steo  where  TRAT  is  converted  to  RDX 
requires  a  slightly  stronger  nitration  reagent  (HN03/P205)  than  the  HNQ3/ 

polyphosphoric  acid  used  in  the  HMX  synthesis  to  give  RDX.  The  three  step 
synthesis  of  RDX  is  diagramed  in  Figure  1 1  and  described  in  detail  in  th9  text 
below. 


1 .3.5-triacety l-hexahydro- 1 .3 .5-triazi ne  (TRAT).62-64  Hexamethylene¬ 
tetramine,  8.3  g  (0.059  mol),  was  added  to  34.3  g  (0.34  mcl)  of  acetic  anhydride 
at  room  temperature  with  stirring.  The  resulting  solution  was  heated  at  98  °C  for 
two  hours.  This  being  complete,  the  reaction  mixture  was  cooled  to  5  °C,  20  mL 
of  distilled  water  was  added,  and  the  resulting  mixture  stirred  for  30  minutes. 
The  solution  was  then  reduced  to  a  viscous  yellow  liquid  under  vacuum. 
Distilled  water  (20  mL)  was  added,  the  mixture  was  cooled,  and  seed  crystals 
were  introduced  to  induce  precipitation.  The  solid  product  was  then  filtered  and 
dried  in  vacuo  over  NaOH  pellets. 

Yield  8.65  g  (64%),  m.p.  95-96  °C,  ’H  NMR  (DMSO-d6)  5  2.1  (-CHa),  5 
5.2  (-CHZ-). 
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FIGURE  11  RDX  Synthesis. 


1.3.5-triac.9tyl-hflxahydrfl,:l>a.S-.lfiaziQa:Cfe-lXBAJ-Jfe)  64  The  trat-c/6 

was  synthesized  using  7.00  g  (0.046  mol)  ot  HMTA-d12  and  28.7  g  (0.28  mol)  of 

acetic  anhydride  in  the  same  manner  described  for  the  TRAT  synthesis  above. 
Yield  7.60  g  (75%),  m.p.  95-96  °C,  NMR  of  TRAT-d6  showed  greater 

than  98%  deuterium  substitution  on  the  ring  methylene  groups  calculated  from 
the  acetyl  methyl  protons. 

Hexahvdro-1 .3.5-trinitro-s-triazine  (RDXl5962'64  As  with  the  final  nitration 
step  in  the  synthesis  of  HMX,  the  composition  of  the  atmosphere  above  the 
nitrating  solution  significantly  affects  product  yield.  An  apparatus  identical  to  the 
one  used  for  the  nitration  of  DADN  was  employed  for  the  nitration  of  TRAT. 

To  5.4  g  (0.378  mol)  of  phosphorus  pentoxide  (weighed  out  under 
nitrogen)  in  a  three-neck  5C0-mL  round-bottom  flask  was  added  dropwise  108 
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mL  of  anhydrous  nitric  add.  The  addition  of  the  HN03  caused  an  exotherm  to 

41  °C.  An  overhead  mechanical  stirrer  was  utilized  to  ensure  complete  mixing 
of  the  HN03  and  P2Os.  The  mixture  was  cooled  to  36  °C  and  5.4  g  (0.026  mo!) 

of  TRAT  was  added  rapidly.  The  reaction  flask  was  placed  In  an  oil  batn  and 
heated  at  65  °C  for  15  minutes.  As  the  temperature  increased,  the  atmosphere 
above  the  reaction  mixture  changed  from  a  light  yellow  to  a  dark  red-brown. 
The  reaction  mixture  was  allowed  to  cool  to  room  temperature  and  then  poured 
on  ice.  After  the  ice  melted  the  white  precipitate  was  filtered,  washed  with 
water,  and  allowed  to  air  dry. 

Yield  4.00  g  (71%).  m.p.  202-203  °C.  ’H  NMR  (DMSO-cfe)  5  6.1  (-CH2-). 

Hexahydro-1JL5-trmilro-s-triazine-dfl  (RDX-dflt.84  RDX-de  was 
synthesized,  using  5.55  g  (0.0253  mol)  of  TRAT-cf6,  following  the  procedure 
described  for  the  synthesis  of  RDX. 

Yield  4.00  g  (69%),  m.p.  203-204  °C,  ’H  NMR  spectrum  showed  greater 
than  98%  isotopic  integrity. 

HMX  and  RDX  Particle  Size 

The  90%.  50%,  and  10%  particle  diameters  and  the  calculated  specific 
surface  area  for  each  of  the  four  synthesized  nitramines  are  listed  in  Table  IV.65 
The  particle  diameter  data  were  obtained  using  a  Malvern  3600E  Particle  Size 
Analyzer.  The  Malvern  analyze  collects  data  by  using  a  light  scattering 
technique  to  determine  the  number  of  particles  in  each  of  a  large  number  of 
diameter  bands.  Using  these  data  the  diameters  listed  in  Table  IV  were 
calculated.  An  example  of  how  the  data  in  Table  IV  are  interpreted  is  as 
follows:  the  particle  size  listed  for  the  90%  diameter  is  simply  the  diameter 
below  wh->h  90%  of  the  sample  particle  diameters  fall.  The  Malvern  determines 
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the  surface  area  of  each  band  using  the  median  diameter  of  the  band  and 
assuming  that  all  particles  are  perfect  spheres.  The  specific  surface  area  for  a 
sample  is  simply  the  sum  of  the  surface  areas  for  all  the  diameter  bands. 
Because  many  particle  samples  do  not  possess  a  normal  distribution  of  sizes,  it 
is  perfectly  conceivable  that  two  samples  with  different  particle  distributions  can 
have  the  same  specific  surface  area. 


TABLE  IV 

Particle  Diameters  (microns,  pm)  and 
Calculated  Specific  Surface  Area  (m2/cm3) 


Compound 

90% 

50% 

to% 

Specific  Surface  Area 

HMX 

47.4 

21.9 

6.2 

0.13 

HMX-da 

60.3 

35.7 

8.6 

0.10 

RDX 

61.2 

36.7 

11.8 

0.10 

RDX-dg 

53.7 

22.8 

8.3 

0.11 

DSC  EXPERIMENTS 

The  p  -►  8  phase  transition  of  HMX  and  HMX-dg  was  investigated  using  a 
Porkin-Elmer  Model  DSC-4  Differential  Scanning  Calorimeter.  Experiments 
consisted  of  placing  the  sample  (either  5.0  or  10.0  ±  0.2  mg)  in  a  sealed 
aluminum  sample  pan,  Perkin-Elmer  Part  No.  219*0041,  and  then  performing  a 
normal  DSC  experiment  scanning  from  150  to  200  °C  at  5  °C/min.  Following 
the  completion  of  the  experiment,  peak  analysis  was  performed  on  the  phase 
transition  endotherm  to  determine  the  onset  temperature,  the  peak  temperature, 
and  the  heat  of  transition. 
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COMMERCIALLY  OBTAINED  CHEMICALS 


The  Ingredients  used  in  the  model  propellant  formulations  were  standard 
propellant  grade  unless  otherwise  noted  and  were  obtained  from  the  following 
manufacturers:  dibutyttln  diacetate  (DBTDA),  reagent  grade,  Alfa  Products; 
hydroxyl-terminated  polybutadiene  (HTPB,  trade  name  R-4SM)  Lot  #803445, 
ARCO;  isophorone  diisocyanate  (IPDI)  Lot  #6932-0008,  Chemischewerke  HQ  Is 
AG;  poly(diethyleneglycol  adipate)  hydroxyl-terminated  (trade  name  R-18)  Lot 
-#D-018  4-002,  Mobay;  trimethylolethane  trinitrate  (TMETN)  Lot  #6H-10M, 
Trojan;  triphenyl  bismuthine  (TPB),  reagent  grade,  Aldrich. 

For  experiments  not  requiring  isotopically  labelled  nitramines,  propellant 
grade  HMX  and  RDX  were  used.  The  HMX,  Lot  #78L675-014  and  weight 
median  diameter  of  26  pm,  and  the  RDX,  Lot  #77H1 50-058  and  weight  median 
diameter  of  35  pm,  were  obtained  from  Aerojet  Corporation. 

PRESSED  PELLET  PREPARATION 

The  bottom  halves  of  No.  3  gelatin  capsules  (Eli  Lilly  Co.)  were  used  to 
contain  the  pressed  nitramine  powder  and  to  act  as  an  inhibitor  to  ensure  the 
uniform  horizontal  burning  of  the  pellets.  The  pellets  were  prepared  by  first 
placing  the  gelatin  capsule  into  the  pellet  die  (Figs.  12  and  13)  then  filling  the 
capsule  with  the  appropriate  nitramine  powder.  The  powder  was  then  lightly 
compressed  to  allow  for  the  addition  of  more  nitramine;  more  powder  was 
added  and  then  compacted.  This  compacting  procedure  was  repeated  several 
times  until  the  level  of  the  pressed  powder  was  a  few  millimeters  above  the  top 
lip  of  the  capsule.  A  hydraulic  press  (Model-C,  Fred  S.  Carver  Inc.)  and  a  punch 
(dia.  0.208  in),  which  fit  snugly  Inside  the  gelatin  capsule  (I.D.  0.211  in),  was 
used  to  compress  the  powder  to  a  pressure  of  50  psi  on  a  1-7/8  inch  ram  (Fig. 
12).  The  punch  was  then  withdrawn  and  more  powder  added  to  a  level  above 
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PELLET  PRESS 


FIGURE  12.  Pellet  Pressing  Apparatus. 
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FIGURE  13.  Pellet  Die,  Holder,  and  Capsule  Support. 
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the  capsule  lip.  The  powder  was  then  pressed  to  its  final  pressure  of  100  psi  on 
the  1-7/6  inch  ram.  This  final  pressure  equates  to  a  pressure  on  the  0.211 -inch 
diameter  pellet  of  7.900  psi.  The  powders  were  pressed  to  this  pressure  in 
order  to  give  pellet  densities  approximately  the  same  as  the  densities  of  the 
formulated  propellants. 

The  pellet  was  removed  by  placing  the  die  with  capsule  onto  the  die 
holder  (Fig.  13)  and  pressing  the  capsule  out  the  bottom  of  the  die  with  the 
punch.  A  razor  blade  was  used  to  shave  off  the  excess  powder  down  to  the  top 
lip  of  the  gelatin  casing.  The  densities  of  the  RDX,  RDX-cfe,  HMX,  HMX-d8 

pellets  are  listed  in  Table  VIII  along  with  the  densities  of  the  propellant  pellets. 
After  pressing,  the  pellets  were  approximately  13  mm  by  5.4  mm. 


MODEL  PROPELLANT  PREPARATION 


Propellant  Formulations 

The  nitramine  propellant  formulations  used  in  this  study,  designated 
CW5,  PB,  and  RRHR-18,  are  listed  in  Tables  V,  VI,  VII  respectively.  The  CW5 
formulation  (Table  V)  was  chosen  because  of  its  simplicity  and  because  it 
contained  both  an  energetic  plasticizer  and  a  high  oxygen  content  polymer.66 
The  CW5  formulation  is  a  relatively  energetic  nitramine  propellant  which  burns 
very  uniformly  owing  to  an  almost  equal  balance  between  fuel  and  oxidizer 
components. 

The  PB  propellant  formulation  was  also  chosen  for  its  simplicity.  In 
ad  lition,  it  has  been  used  previously  in  studies  involving  both  combustion 
modeling67  and  propellant  burning  properties.68  Thus,  a  large  amount  of 
information  exists  dealing  with  the  PB  propellant  burning  characteristics.  The 
PB  propellant  is  much  less  energetic  than  the  CW5  propellant.  It  has  a  fuel  rich 
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formulation  which  bums  in  a  very  heterogeneous  manner  giving  off  a  great  deal 
of  carbonaceous  residue.  The  CW5  and  PB  propellants  represent  the  opposite 
extremes  as  far  as  energetics  and  burning  properties,  and  therefore  they  should 
provide  a  relevant  test  of  isotope  effects  in  cured  nitramine  propellants. 


TABLE  V 

CWS  Propellant  Formulation 


Comoound 

Percent 

Nitramine 

73.20 

R-18  (polymer) 

0  0 

II  li 

8.51 

HOHC^CHgOCtCH^COifl-H 

TMETN  (plasticizer) 

ch2ono2 

HgC-C— CH20N02 

ch2ono2 

16.89 

IPDI  (curative) 

H,C  _  NCO 

1.39 

DBTDA  (cure  catalyst) 

Bu2Sn(OAc)2 

0.01 

The  RRHR-18  propellant  formulation  was  not  used  in  any  isotope  effect 
experiments.  It  was  designed  to  provide  an  additional  propellant  to  study  the 
effects  of  combustion  pressure  and  propellant  formulation  on  flame  structure 
and  burning  surface  characteristics. 
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TABLE  VI 


Percent 

Nitramine 

80.00 

R-45M  (polymer) 

18.52 

HO — |(CHtCH=sCHCHJ^-(CH,CH,) — (CH,CH=CHCH,>- 

OH 

L  CH  =CH, 

n 

IPDI  (curative) 

1.46 

TPB  (cure  catalyst) 

0.02 

(C6H5)3Bi 

TABLE  VII 

ComDOund 

Percent 

Nitramine 

80.00 

R-18  (polymer) 

17.17 

IPDI  (curative) 

2.81 

DBTDA  (cure  catalyst) 

0.02 

Propellant  Binders 

In  our  experiments  two  different  types  of  binders  were  used.  The  binder 
of  the  CW5  ano  RRHR-18  formulations  employed  a  polyester  polymer  R-18 
which  was  37%  by  weight  oxygen  and  had  a  number  average  molecular  weight 
of  1400.  The  CW5  formulation  also  employed  a  nitrate  ester  plasticizer 


,  ,h.  snef0V  ot  the  binder  in  addition  to  Increasing 

^  Th.  PB  formulation  used  a  pol, butadiene  polymer  B-45M 

binder  elasticity.  averafle  molecular  weight  ot 

wtilch  was  3.6%  by  weight  oxygen  wrth  nifl«rticizer 

20qo  Because  th,  PB  and  RRHR*1 8  formulations  did  no,  contain  a  plasticiz  . 

„*r  binder  wa.  much  harder  than  thaf  in  ,h.  CW5  —on. 


..  wm 


H$r  CHjNCO 

IPDI 


R-45M 


alcohol 

ISOCYANATE 

COMPLEX 


FIGURE  14.  urethane  Cure  of  the  R-45M  Polymer. 
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Both  polymers  were  crosslinked  using  a  so  called  "urethane  cure."  This 
involves  the  condensation  reaction  ot  the  isocyanate  curative  with  the  hydroxyl 
end  group  on  the  polymer  to  form  a  urethane  linkage  (Fig.  14).69  The  cure 
catalyst  appears  to  facilitate  the  curing  process  by  stabilizing  the  alcohol- 
isocyanate  complex  thus  enhancing  the  rate  of  formation  of  the  urethane 
linkages.70 

The  elastomeric  binder  consists  of  a  three  dimensional  matrix  of  polymer 
units  linked  together  by  the  isocyanate  curative.  The  amount  of  curative  used 
dictates  the  amount  of  crosslinking  in  the  binder  and  thus  significantly  affects 
the  bulk  properties  of  the  binder.  In  the  caso  of  the  CW5,  RRHR-10,  and  PB 
formulations,  a  one-to-one  mole  ratio  of  polymer  to  curative  was  used. 

CAUTION:  Care  must  be  taken  while  mixing  and  handling  the  HMX  and  RDX 
propellants  due  to  their  explosive  nature.  Solid  Teflon  utensils  and  beakers 
were  used  at  all  times  in  the  propellant  preparation  to  preclude  accidental  spark 
ignition. 

Propellant  Mixing 

The  propellants  were  prepared  on  a  2.5-g  scale  in  a  50-mL  Teflon 
beaker.  Initially,  the  polymer,  curative,  and,  if  necessary,  plasticizer  w  e  mixed 
into  a  viscous  homogeneous  liquid.  The  nitramine,  which  had  been  dried  at  60 
°C  for  12  h,  was  then  mixed  in  until  a  smooth,  pasty  consistency  was  achieved. 
Finally,  the  cure  catalyst  was  added  and  the  completed  propellant  carefully 
packed  into  the  bottom  halves  of  five  No.  3  gelatin  capsules.  Great  care  was 
taken  to  avoid  leaving  any  trapped  air  in  the  propellant.  The  propellant-filled 
capsules  were  then  cured  at  40  °C  in  a  pressure  bomb  (Parr  Model  1108 
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Oxygtn  Bomb)  pressurized  to  50  atm  with  nitrogen.  The  CW5  and  RRHR-18 
propellants  were  fully  cured  after  12  h  at  40  °C.  However,  the  PB  propellants 
required  from  24  to  36  h  at  40  °C  to  fully  cure.  The  density  of  the  propellants  are 
listed  in  Table  VIII.  The  cured  CW5,  RRHR-18,  and  PB  propellant  pellets  were 
approximately  13.5  mm  by  5.4  mm. 


TABLE  VIII 


HMX  /  HMX-Cfe 

1.525  ±0.015 

RDX/RD  X-flfe 

1.400  10.010 

HMX  /  HMX-c^  CW5 

1.616  ±0.015 

RDX  /  RDX-d^  CW5 

1.480  1  0.030 

HMX  /  HMX-cfe  PB 

1.507  1  0.016 

RDX/RDX-cfe  PB 

1.39510.025 

HMX  RRHR-18 

1.554  1  0.0* ' 

RDX  RRHR-18 

1.42310.0.0 

The  reported  error  is  at  the  95%  confidence  level. 


PVC  Coated  Propellants 

The  gelatin  capsules  frequently  cracked  due  to  dehydration  when 
propellant  curing  required  more  than  24  h.  Because  the  propellant  burning 
surface  propagated  down  these  cracks,  making  accurate  burn  rate 
determination  almost  impossible,  the  damaged  gelatin  capsules  did  not  function 
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effectively  as  burning  inhibitors.  Therefore,  several  experiments  were 
performed  using  a  polyvinyl  chloride  (PVC)  coating  as  the  burning  inhibitor. 
This  was  accomplished  by  first  preparing  the  propellant  in  the  gelatin  capsules 
as  described  above.  The  fully  cured  encapsulated  propellants  were  then 
placed  in  a  dilute  aqueous  acetic  acid  solution  to  dissolve  away  the  gelatin. 
The  propellant  pellets  were  then  thoroughly  washed  with  distilled  water  and  air 
dried.  The  pellets  were  dipped  into  a  dilute  solution  of  PVC  in  methylene 
chloride  and  air  dried.  Finally,  a  razor  blade  was  used  to  shave  off  a  small  disk 
of  propellant  from  the  top  of  the  pellet  to  give  an  uncoated  surface  for  burning. 

WINDOW  BOMB  COMBUSTION  PROCEDURE 

Burn  rate  experiments  were  performed  in  a  reinforced  concrete  test  cell 
at  the  Air  Force  Astronautics  Laboratory,  Propellant  Experimental  Area  1-30.  All 
experiments  were  monitored  remotely  from  a  control  room  adjacent  to  the  test 
cell. 

Window  Bomb 

The  window  bomb  (AFRPL  Window  Bomb  drawing  #X7616176)  used  for 
all  combustion  experiments  is  shown  in  Figure  15.  Samples  in  the  window 
bomb  were  illuminated  using  a  xenon  spot  light  (Color  Arc  2000).  Experiments 
were  recorded  using  a  high  speed  video  system  (Ektapro  1000,  Kodak-Spin 
Physics).  Data  was  collected  at  a  rate  of  10C0  frames  per  second  with  timing 
data  recorded  on  each  video  frame. 

Burn  Rate  Experiments 

A  typical  burn  rate  experiment  began  by  first  placing  a  ruler  on  the  base 
of  the  window  bomb  in  the  same  location  as  occupied  by  the  pellet  (Fig  16a). 
The  base  was  then  screwed  on  to  the  bomb  and  approximately  one  hundred 
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FIGURE  15.  High  Pressure  Window  Bomb. 
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BASE  OF  THE 
WINDOW  BOMB 


FIGURE  16,  Placement  of  Pellet  and  Scale. 


frames  of  the  scale  were  taken  using  the  high  speed  video  system.  These  data 
were  later  used  to  determine  burn  rates.  The  scale  was  then  removed  and 
replaced  with  a  pellet  placed  atop  a  small  piece  of  putty  (Fig.  16b).  The  pellet, 
with  a  Nichrome  wire  stretched  across  the  top,  and  the  base  were  then  attached 
to  the  bomb.  Normally  the  window  bomb  was  pressurized  to  1000  psi  with 
nitrogen  prior  to  the  experiment.  The  pressure  inside  the  window  bomb  was 
regulated  to  ±  10  psi  during  the  experiment.  A  constant  15  psid  (a  15  psi 
differential  between  inlet  and  exit  pressures)  nitrogen  purge  during  the 
experiment  maintained  a  clear  photographic  path  to  the  pellet  by  removing 
combustion  products  and  dissipating  convective  heat  currents.  The  pellet  was 
ignited  by  the  Nichrome  wire  and  the  burning  recorded  using  the  high  speed 
video  system. 

Bum  rates  were  determined  in  the  following  manner.  A  single  video 
frame  of  the  ruler  inside  the  window  bomb  was  displayed  on  a  video  monitor.  A 
clear  mylar  sheet  containing  a  grid  which  exactly  matched  the  ruler  markings 
displayed  on  the  monitor  was  taped  to  the  screen.  The  video  of  the  pellet  burn 
was  then  advanced  until  all  the  effects  of  ignition  had  dissipated.  This  was 
normally  after  about  the  top  20%  of  the  pellet  had  burned.  At  this  point  the  time 
given  on  the  video  frame  and  the  distance  on  the  grid  were  rocorded.  The 
video  was  then  advanced  until  the  pellet  burn  reached  the  point  where  the 
curvature  of  the  pellet  bottom  began.  At  this  location  the  time  and  distance  were 
again  recorded.  The  pellet  burn  rate  was  calculated  simply  by  taking  the 
distance  burned  over  the  change  in  time. 

Extinguishing  Burning  Propellants 

In  order  to  observe  the  composition  and  structure  of  the  propellant 
burning  surface,  several  burning  pellets  were  extinguished  through  rapid 
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decompression  of  the  combustion  chamber.  To  begin,  the  procedure  followed 
exactly  that  used  for  normal  pellet  bums.  However,  once  the  pellet  was  ignited, 
the  combustion  chamber  was  vented  to  the  outside  atmosphere,  causing  rapid 
decompression  and  extinguishing  the  burn.  Decompression  from  1000  psi  to 
ambient  pressure  required  approximately  one  second. 

Electron  micrographs  of  the  propellant  burning  surfaces  were  obtained 
using  an  ISI  CL-6  scanning  electron  microscope  (SEM).  Samples  were  coated 
with  a  gold  film  prior  to  SEM  examination  using  a  Hummer  X  sputter  coater  to 
facilitate  the  efficient  conduction  of  electrons  away  from  the  propellant  surface. 
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RESULTS  AND  DISCUSSION 


The  linear  burn  rate  of  nitramines  and  nitramine  propellants  depends  on 
both  kinetic  and  physical  parameters. 3,88  In  our  experiments,  factors  such  as 
pressure,  temperature,  the  shape  and  size  of  the  burning  surface,  etc.  are 
essentially  the  same  for  two  samples  differing  only  by  isotopic  substitution. 
Thus  any  burn  rate  differences  should  be  due  to  the  effect  of  isotopic 
substitution  on  the  rate  limiting  process  in  the  combustion  mechanism. 
Because  it  is  virtually  impossible  to  obtain  kinetic  rate  constants  from 
combustion  experiments,  burn  rates  must  be  used  to  obtain  isotope  effect 
values.  Thus,  to  make  any  conclusions  about  the  rate  limiting  process  from 
KDIE  values  one  must  understand  the  relationship  between  burn  rate  and  the 
kinetic  rate  constant  of  the  combustion  reaction. 

Nitramine  and  nitramine  propellant  combustion  begins  with  condensed 
phase  decomposition  reactions  at  the  burning  surface  and  proceeds  through  to 
gas  phase  oxidation  reactions  in  the  luminous  flame.3,18,71  Reaction 
temperatures  for  combustion  range  from  approximately  673  K  at  the  burning 
surface  to  in  excess  of  1800  K  in  the  luminous  flame.68  A  fundamental  principle 
of  reaction  kinetics  is  that  reaction  rates  increase  with  increasing  temperature/3 
In  addition,  it  can  be  shown  that  reaction  rates  tend  to  increase  when  going  from 
a  non-mobile  condensed  phase  to  the  gas  phase.74  Thus,  it  would  seem,  a 
priori,  that  the  slowest  reactions  in  the  combustion  mechanism  occur  in  the 
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condensed  phase.  The  numerous  mechanistic  models  proposed  for  HMX  and 
RDX  condensed-phase  decomposition  suggest  a  limited  number  of  reactions 
occurring  in  the  condensed  phase  due  to  the  early  production  of  gas-phase 
intermediates.8*8*11-19-75  Recent  theoretical  studies  indicate  that  after  the  early 
bond  breaking  steps  of  HMX  and  RDX  decomposition,  the  further  breakup  of  the 
nitramine  becomes  relatively  facile.78  From  this,  one  can  reasonably  conclude 
that  it  is  likely  that  the  slowest  step  in  the  nitramine  combustion  mechanism 
occurs  early  in  the  condensed-phase  nitramine  decomposition. 

Depending  on  the  initial  conditions,  the  decomposition  of  HMX  and  RDX 
follow  either  a  first  order  (eq.  24)17*20*77*78  or  an  autocatalytic  rate  law  (eq.  25). 19 

k(1  -  a)  (24) 

~-ka(1  -  a)  (25) 

The  extent  of  reaction  is  represented  by  a,  t  is  time,  and  k  is  the  kinetic  rate 
constant  for  the  reaction.  It  can  be  shown  that 

— -=rC  (26) 

where  r  is  the  linear  burn  rate  and  C  is  a  shape  factor  which  depends  on  the 
dimensions  of  the  burning  surface.79  Assuming  the  same  process  that  controls 
the  rate  of  decomposition  controls  the  rate  of  combustion,  then  with  respect  to 
the  equations  above  (eq.  24,  25,  and  26),  the  linear  burn  rate  should  be 
proportional  to  the  kinetic  rate  constant  such  that  KDIE  values  can  be  obtained 
from  burn  rates  by  using  the  following  equation. 


We  have  assumed  equation  24  to  be  valid,  and  have  used  it  to  calculate  all 
KDIEs  presented  in  this  tliesis. 

The  rate  constant  k  for  the  combustion  reaction  is  dependent  on  the  rate 
constant  of  the  slowest  step  and  to  a  lesser  extent  on  the  steps  preceding  the 
slow  step.73  Since  the  rate  controlling  process  for  nitramine  combustion 
appears  to  occur  early  in  the  decomposition,  the  value  of  k  should  reflect 
primarily  the  rate  <  »i.wt  of  the  slow  step.  Thus,  the  KDIE  values  obtained 
from  burn  rates  should  be  reasonably  representative  of  the  effect  of  deuterium 
substitution  on  the  rate  controlling  process 

Results  from  the  computational  modeling  of  nitramine  propellant  burn 
rates  suggest  that  k  is  proportional  to  r2  instead  of  r.87  00'85  However,  these 
combustion  modeling  studies  are  generally  empirical  in  nature  and  only  loosely 
based  on  kinetic  data.  Therefore,  it  is  difficult  to  attach  any  significance  to  the 
implications  of  the  modeling  results.  Even  if  k  were  proportional  to  r2.  this  would 
simply  mean  that  the  true  KDIE  values  for  the  nitramine  and  nitramine 
propellant  combustion  reaction  rate  constants  would  be  some  what  larger  than 
the  values  shown  in  this  thesis  (in  particular,  the  square  of  the  KDIEs  based  on 
r).  In  the  end  this  would  not  significantly  affoct  the  interpretations  we  have 
made. 

BURN  RATE  EXPERIMENTS 

It  was  the  intent  of  the  following  experiments  to  measure  the  relative  burn 
rate  differences  between  samples  containing  either  the  protiated  or  deuteriated 
nitramine.  Additionally,  it  was  our  purpose  to  compare  the  KDIEs  of  the  pure 
nurammes  with  the  KDIEs  of  the  corresponding  propellant  mixtures  Toward 
thi  i  end,  every  effort  was  made  to  preclude  differences  in  burn  rate  due  to 
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synthesis,  sample  history,  propellant  mixing,  pellet  preparation,  and  pellet 
density. 


Pressed  Pellet  Results 

As  a  baseline  for  future  comparison  to  model  nitramine  propellants  we 
duplicated  the  previous16*17  combustion  studies  on  HMX,  HMX-d8,  RDX,  and 
RDX-d6  pressed  pellets,  acquiring  the  burn  rate  data  listed  in  Table  IX.  The 
burn  rates  of  commercial  grade  HMX  and  RDX  pressed  pellets  were  also 
measured  giving  values  of  0.535  ±  0.057  in/s  and  0.551  1  0.020  in/s, 
respectively.  These  burn  rates  are  in  the  same  range  as  the  values  listed  lor  the 
synthesized  HMX  and  RDX  indicating  the  results  given  in  Table  IX  are 
reasonable. 


TABLE  IX 

Bum  Rate  Results  for  HMX,  HMX-dg,  RDX,  and 
RDX-d8  Pressed  Pellets  at  100C  psi." 


HMX  <rH)t 

HMX-fll  (r0)» 

RDX  (rH)» 

RDX-<%(rD)» 

0.503 

0.394 

0.600 

0.449 

0.505 

0.386 

0.587 

0.456 

0.531 

0.395 

0.609 

0.457 

0.503 

0.397 

0.578 

0.452 

rH  -0.51010.022 

r0  -  0.393  ±  0.008 

rH  -  0.594 1 0.022 

r0  -  0.454  1  0.006 

(}  y 

KDIE-  f- 

\  ») 

-  1.3010.04 

(t  } 

KDIE-  J1 

It V 

»  1.31  10.03 

•  Reported  error  is  at  the  95%  confidence  level.  b  Values  reported  in  in/s. 
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The  KDIE  of  1 .40  ±  0.1 1  for  HMX  combustion  at  varying  pressures  (500- 
1500  psi)  reported  previously  by  Shackelford  et  al.16  is  consistent  with  our  95% 
confidence  interval  of  1 .30  ±  0.04  for  the  mean  isotope  effect  for  HMX.  The 
previous  report  of  1.46  ±  0.23  for  RDX17  under  conditions  identical  to  ours  (1000 
psi)  is  also  consistent  with  our  estimate  of  1 .31 1 0.03  for  pure  RDX. 

Using  the  evaluation  criteria  developed  in  the  previous  chapter,  the 
estimate  for  the  isotope  effect  arising  from  a  C-H  bond  rupture  in  HMX  or  RDX  is 
1 .41  (see  page  37).  The  estimated  isotope  effect  for  a  C-N  bond  rupture  is  i  .19 
(see  page  37).  These  estimates  were  made  from  data  collected  at  temperatures 
around  500  K,  whereas  the  temperature  of  the  burning  surface  during 
combustion  is  approximately  673  K.  If  one  used  equation  16  and  assumed  that 
the  isotope  effects  were  due  strictly  to  zero  point  energy  differences  (i.e.,  (aE8]q), 

then  a  temperature  extrapolation  can  be  performed  on  the  above  estimates 
giving  KDIE  values  of  1.29  for  a  C-H  bond  rupture  and  1.14  for  a  C-N  bond 
rupture  at  673  K. 

The  rate  limiting  process  for  the  decomposition  of  HMX  and  RDX  under 
most  circumstances  appears  to  be  a  C-H  bond  rupture.1 7,19,20  Assuming  a 
consistent  rate  limiting  process  for  both  decomposition  and  combustion  and 
applying  the  above  KDIE  criteria  for  673  K,  the  isotope  effect  values  for  both 
HMX  and  RDX  of  1.30  ±  0.04  and  1.31  ±  0.03,  respectively  indicate  a  rate 
limiting  OH  bond  rupture  in  the  combustion  process.  These  assignments  are  in 
agreement  with  the  previous  HMX  and  RDX  combustion  studies.10,17 

Effects  of  Curing  Conditions  and  Burning  Inhibitors 

The  effects  of  curing  conditions  and  burning  inhibitors  on  propellant  burn 
rate  were  investigated  in  order  to  test  the  ruggedness  of  our  CW5  propellant 
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formulation.  The  results  of  these  experiments  on  CW5  propellant  formulations 
containing  commercial  grade  HMX  or  RDX  are  listed  in  Table  X. 


TABLE  X 


The  Effects  of  Burning  Inhibitor,  Curing,  and  Curing  Pressure 
on  CW5  Propellant  Burn  Rates  at  1000  psi. 


Formulation 

Cure  Conditions8 

Burning  Inhibitor 

Bum  Rate  (in/s)1 

HMX-CW5 

1  atm 

Gelatin  Capsule 

0.175  ±0.006 

HMX-CW5 

50  atm 

Gelatin  Capsule 

0.166  ±0.003 

HMX-CW5 

uncuredc 

Gelatin  Capsule 

0.158  ±0.007 

HMX-CW5 

50  atm 

PVC 

0.177  ±0.004 

RDX-CW5 

1  atm 

Gelatin  Capsule 

0.177  ±0.003 

RDX-CW5 

50  atm 

Gelatin  Capsule 

0.178  ±0.005 

RDX-CW5 

uncured® 

Gelatin  Capsule 

0.174  ±0.004 

RDX-CW5 

50  atm 

PVC 

0.186  ±  0.009 

*  All  propellants  were  maintained  at  40  °C  until  fully  cured.  b  Reported 
error  is  at  the  95%  confidence  interval.  c  No  DBTOA  cure  catalyst  was 
added  to  these  propellants. 


Large  scale  propellant  mixing  operations  avoid  the  problem  of  trapped 
gases  by  performing  most  propellant  mixing  procedures  under  a  vacuum. 
However,  because  it  is  almost  impossible  to  use  normal  mixing  methods  when 
working  on  a  2.  5  g  scale,  our  initial  propellant  mixtures  were  prepared  and 
cured  at  atmospheric  pressure.  Not  unexpectedly,  these  propellants  contained 
trapped  gases  in  the  form  of  bubbles  in  the  propellant.  When  pellets  ot  those 
propellants  wore  placed  inside  the  window  bomb  and  then  pressurized  to  1000 


77 


psi  the  sides  of  the  pellets  were  compressed  due  to  the  compression  of  the 
bubbles  inside  the  cured  propellant.  This  problem  was  avoided  if  the  propellant 
was  cured  at  50  atm  (735  psi).  Evidently,  the  increased  curing  pressure  either 
forced  out  or  compressed  the  bubbles  in  the  propellant.  As  one  can  see  from 
the  data  in  Table  X  the  effect  of  the  change  in  curing  pressure  is  very  small. 
There  was  only  a  slight  decrease  in  the  HMX-CW5  bum  rate  and  no  change  in 
the  RDX-CW5  bum  rate  when  changing  the  cure  pressure  from  1  to  50  atm. 

The  difference  in  the  bum  rate  of  cured  versus  uncured  HMX-CW5  is 
relatively  small.  In  the  case  of  the  RDX-CW5  propellant  the  difference  is 
insignificant.  These  results  indicate  the  relative  unimportance  of  the  physical 
structure  of  the  propellant  binder  in  determining  burn  rate.  Consequently,  any 
differences  in  crosslinking  density,  and  thus  binder  physical  properties, 
between  identical  propellant  mixes  should  not  introduce  any  significant 
variation  in  bum  rate. 

There  was  a  small  increase  in  the  bum  rate  of  both  HMX  and  RDX  CW5 
propellants  when  the  burning  inhibitor  was  switched  from  a  gelatin  capsule  to 
PVC.  More  than  likely  these  differences  are  due  to  the  different  physical 
characteristics  of  the  inhibitors.  The  gelatin  capsule  was  better  at  limiting  the 
propellant  burning  to  the  top  surface  of  the  pellet.  The  PVC  coating  allowed  the 
propellent  burning  to  proceed  down  the  side  of  the  pellet  slightly  ahead  of  the 
advancing  top  surface;  this  allowed  more  burning  surface,  thus  increasing  the 
overall  linear  burn  rate. 

Overall  one  can  see  from  the  data  in  Table  X  that  the  variation  in 
parameters  had  only  a  minor  effect  on  the  CW5  propellant  bum  rate.  Thus  one 
can  conclude  that  any  changes  in  burn  rate  between  normal  and  deuterium 
labeled  propellants  should  be  due  strictly  to  the  effect  of  deuterium  on  the  rate 
controlling  process. 
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TABLE  XI 


Bum  Rate  Results  for  HMX  and  HMX-d8  CW5  Propellants  at  1 000  psi.a 


HMX-CW5  (rH)b 

HMX-cfe*CW5  (r0)b 

0.173 

0.180 

0.142 

0.174 

0.179 

0.141 

0.178 

0.201 

0.131 

0.159 

0.155 

0.137 

0.162 

0.206 

0.136 

rH- 0.177  ±0.012 

fD- 0.137  ±0.005 

f  O 

KDIE-  H 

VI Dj 

-1.29  ±0.07 

•  Reported  error  is  at  the  95%  confidence  level.  ^  Values  reported 
in  in/s. 


KDIE  in  CW5  Propellant  Mixtures 

The  burn  rate  results  for  normal  and  deuterium  labeled  HMX  and  RDX 
CW5  propellants  are  listed  in  Tables  XI  and  XII  respectively.  Because  of  the 
small  size  of  the  KDIE  for  data  Set  I  and  the  low  burn  rate  of  the  RDX-CW5 
propellant  when  compared  with  the  burn  rate  of  the  commercial  grade  RDX* 
CW5  (Table  IX),  a  second  set  of  data  was  obtained.  The  burn  rates  for  both  the 
RDX-CW5  and  RDX-de-CW5  propellants  In  the  second  data  set  were 

significantly  greater  than  in  the  first  data  set.  The  propellants  for  the  second 
data  set  were  mixed  on  a  day  of  unusually  high  humidity  (for  the  desert).  The 
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presence  of  moisture  evidently  caused  the  propellant  to  degas  more  than  usual, 
and  this  apparently  caused  the  higher  overall  bum  rates. 


TABLE  XII 

Bum  Rate  Results  for  RDX  and  RDX-de  CW5  Propellants  at  1000  psi* 


Set! 

Set  II 

RDX-CW5  (rH)b 

RDX-d8-CW5  (r0)b 

RDX-CW5  (rH)b 

RDX-d^-CW5  (rD)b 

0.162 

0.125 

0.241  0.253 

0.193  0.179 

0.164 

0.135 

0.229  0.218 

0.196  0.188 

0.143 

0.135 

0.252  0.245 

0.191  0.188 

0.156 

0.136 

0.227  0.238 

0.173  0.184 

0.158 

0.140 

0.257 

fH«  0.156  ±0.010 

*r0-  0.1 34  ±0.005 

rH- 0.240  ±0.010 

0.1 87  ±0.006 

(f  'i 

KDIE  ■  y- 

W 

-  I.i 6  ±0.05 

(f  \ 
KDIE-  y- 

VO) 

-1.29  ±0.05 

Average  KDIEC 

rn  1 .24  ±  0.05 

a  Reported  error  is  at  the  95%  confidence  level  b  Values  reported  in  in/s. 
c  Weighted  average  for  data  sets  I  and  II. 


It  is  important  to  note  that  each  data  set  is  self  consistent.  In  each  data 
set  both  the  RDX-CW5  and  POX-ds-CW5  propellants  were  prepared  under 

identical  conditions  using  the  same  ingredients,  the  cured  pellets  were  burned 
on  the  same  day  under  identical  conditions,  and  the  data  were  reduced  in  the 
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same  manner.  Therefore,  any  differences  in  bum  rates  between  the  RDX-CW5 
and  RDX-de-CW5  propellants  are  due  strictly  to  the  effect  of  deuterium 

substitution  on  the  rate  limiting  process.  Consequently,  the  differences  in  the 
KDIE  values  for  the  two  data  sets  represent  simple  statistical  variations.  Thus,  it 
should  be  valid  to  average  the  two  KDIE  values. 

The  data  in  Tables  XI  and  XII  clearly  demonstrate  an  isotope  effect  for 
both  HMX  and  RDX  CW5  propellant  mixtures.  The  magnitude  of  the  KDIE  for 
HMX-CW5  of  1.29  ±  0.07  is  almost  identical  to  that  of  the  pure  HMX.  The  KDIE 
for  RDX-CW5  of  1 .24 1 0.05  is  statistically  not  significantly  different  at  the  95% 
confidence  level  from  the  value  for  the  pure  RDX  of  1.31  ±  0  03.  These  data 
provide  evidence  that  the  mechanisms  which  control  the  combustion  of  the  pure 
HMX  and  RDX  continue  to  control  the  combustion  of  the  HMX  and  RDX  CW5 
propellants.  Therefore,  a  C-H  bond  rupture  during  the  condensed  phase 
decomposition  appears  to  be  the  rate  controlling  process  in  the  combustion  of 
both  HMX  and  RDX  CW5  propellant  mixtures. 

KDIE  in  PB  Propellant  Mixtures 

The  burn  rate  results  for  normal  and  deuterium  labeled  HMX  and  RDX 
PB  propellant  mixtures  are  listed  in  Table  XIII.  The  burn  rate  of  PB  propellants 
containing  commercial  grade  HMX  and  RDX  were  measured,  giving  nnan 
values  of  0.1 5^  ±  0.007  In/s  and  0.177  ±  0.006  in/s,  respectively.  The 
commercial  grade  results  appear  to  confirm  the  significantly  lower  HMX-PB 
burn  rate  when  compared  to  the  RDX-PB  and  previous  model  propellant  burn 
rates. 

The  KDIE  of  1.20  ±  0.04  for  the  RDX-PB  propellant  mixture  is  statistically 
different  at  the  95%  confidence  level  from  the  KDIE  value  of  1.31  ±  0.03  for  the 
pure  RDX.  This  decrease  in  isotope  effect  could  be  due  to  a  large  number  of 
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factors.  Unlike  the  CW5  propellant  mixtures,  the  PB  propellants  bum  in  a  very 
heterogeneous  manner  producing  large  amounts  of  carbonaceous  residue. 
This  residue  tends  to  remain  on  the  surface  of  the  burning  propellant,  possibly 
causing  an  increase  in  the  burning  surface  temperature  due  to  increased  heat 
feedback  from  the  flame  to  the  burning  surface.  This  increase  in  the  local 
temperature  of  the  surface  at  which  the  ROX  decomposition  takes  place  would 
be  expected  to  result  in  a  smaller  KDIE  value. 


TABLE  XIII 

Burn  Rate  Results  for  HMX,  HMX-cfe,  RDX  and 
RDX-d^  PB  Propellants  at  1000  psi .■•b 


HMX-PB  (rH)c 

HMX-d8-PB  (rD)c 

RDX-PB  (rH)c 

RDX-de-PB  (r0)'; 

0.124 

0.104 

0.193 

0.169 

0.109 

0.114 

0.189 

0.156 

0.112 

0.110 

0.195 

0.160 

0.116 

0.107 

0.157 

0.111 

0.110 

0.162 

rH«  0.1 14  ±0.007 

r0«  0.109  ±0.005 

rH- 0.192  ±0.008 

f0  «  0.161  ±0.006 

(}  > 

KDIE-  -f 

-1.0510.05 

frj 

KDIE-  — 

-  1.2010.04 

*  Reported  error  is  at  the  95%  confidence  level.  b  PVC  was  used  as  the  burning 

inhibitor  for  all  pellets.  c  Values  reported  in  in/s. 

The  decrease  in  the  magnitude  of  the  KDIE  also  could  be  the  result  of  a  change 
in  the  RDX  decomposition  mechanism.  This  change  might  be  a  new  rate 
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limiting  step  other  than  a  C-H  bond  rupture  (e.g.,  C-N),  or  it  may  be  the 
introduction  of  a  rate  limiting  process  in  addition  to  the  C-H  bond  rupture. 
Finally,  the  difference  in  KOIE  values  may  simply  be  due  to  a  statistical!/ 
inaccurate  estimate  of  the  RDX-PB  propellant  isotope  effect  value  caused  by  a 
small  sample  size. 

The  RDX-PB  KDIE  of  1.20  falls  between  the  primary  and  secondary 
estimates  of  1 .29  and  1.14,  respectively,  at  673  K;  thus  it  is  difficult  to  assign  this 
effect  to  either  a  C-H  or  a  C-N  bond  rupture.  Because  of  the  fact  that  all 
previous  RDX  decomposition  and  combustion  studies  indicate  a  rate  limiting  C- 
H  bond  rupture  and  the  fact  that  the  KDIE  is  larger  than  the  estimated  secondary 
effect,  I  would  venture  to  say  that  the  KDIE  for  the  combustion  of  RDX-PB 
propellants  is  due  to  a  rate  limiting  C-H  bond  rupture.  However,  the  assignment 
is  tenuous. 

The  overall  linear  bum  rates  for  the  HMX-PB  propellant  as  shown  in 
Table  XIII  are  significantly  less  than  for  the  RDX-PB  propellant  and  all  the  CW5 
propellant  mixes.  Furthermore,  the  HMX-PB  propellant  exhibits  a  KDIE  of  only 
1 .05  ±  0.05,  which  is  much  smaller  than  for  all  the  other  propellants  in  this  study. 
These  data  indicate  a  significant  change  in  combustion  mechanism  between 
pure  HMX  and  the  HMX-PB  propellant.  Although  the  isotope  effect  appears  to 
be  real  and  its  value  falls  well  below  the  estimated  secondary  effect  value,  one 
cannot  rule  out  the  possibility  that  the  observed  effect  is  simply  a  very  small 
primary  effect.  Thus,  from  the  data  available  it  is  impossible  to  make  any  kind 
of  assignment  of  a  rate  limiting  process  for  HMX-PB  propellant  combustion. 

Interestingly,  the  "mixed  melt"  region  of  HMX  decomposition  observed  by 
Shackelford  et  al.19  exhibited  an  isotope  effect  very  close  to  unity.  One  cculd 
speculate  that  the  small  KDIE  value  for  the  HMX-PB  propellant  was  due  to  the 
majority  of  HMX  decomposition  occurring  in  the  mixed  melt  as  opposed  to  either 
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the  so Ud  or  liquid  phase.  This  in  turn  could  be  the  result  of  the  heterogeneous 
nature  of  the  PB  propellant  burning. 

PROPELLANT  BURNING  SURFACE  AND  FLAME  CHARACTERISTICS 

Kubota,71  studying  the  combustion  of  a  nitramine  based  propellant 
consisting  of  80%  HMX  or  RDX  and  20%  polyurethane  binder  (28%  by  weight 
oxygen),  found  that  the  burning  propellant  exhibited  a  two  stage  flame 
consisting  of  a  preparation  zone  (dark  zone)  directly  above  the  burning  surface 
and  below  a  luminous  flame  zone.  He  also  found  that  the  thickness  of  this 
preparation  zone  varied  with  pressure  from  approximately  0.3  cm  at  294  psi  (20 
atm)  to  roughly  0.03  cm  at  735  psi  (50  atm). 

In  burn  rate  experiments  with  our  CW5  formulation,  which  differs  from 
Kubota's  propellant  in  that  it  contains  an  energetic  plasticizer  (TMETN)  and  a 
binder  with  a  slightly  higher  oxygen  content  (crosslinked  R-18  polymer,  37%  by 
weight  oxygen),  we  observed  no  preparation  zone  at  pressures  ranging  from 
200  to  1 500  psi.  Figure  1 7  is  indicative  of  the  flame  structure  observed  for  all 
CW5  propellants  over  the  pressure  range  studied. 

A  propellant  (RRHR-18),  identical  to  Kubota's  with  the  exception  that  it 
employed  a  binder  composed  of  crosslinked  R-18  polymer,  was  prepared  and 
burned  at  various  pressures.  There  was  no  significant  difference  between  the 
bum  rates  of  the  RRHR-18  and  CW5  propellants.  However,  at  pressures  below 
750  psi,  the  RRHR-18  propellants  clearly  exhibited  the  two-stage  flame 
observed  by  Kubota.71  Furthermore,  the  preparation  zone  of  the  RRHR-18 
flame  (Fig.  18)  decreased  with  increasing  pressure  in  agreement  with  the 
previous  study.  Additionally,  the  RRHR-18  flame  was  more  diffuse  and 
significantly  lower  in  intensity  than  the  CW5  flame. 

The  position  of  the  luminous  flame  is  a  dynamic  balance  between  the 
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rate  of  production  of  gaseous  reactants  at  the  burning  surface  and  the  rate  of 
progression  of  the  flame  front.72,86  When  the  rate  of  oxidation  reactions  in  the 
luminous  flame  and  thus  the  rate  of  progression  of  the  flame  front  exceeds  the 
rate  of  production  of  gaseous  reactants,  the  luminous  flame  moves  toward  the 
burning  surface  until  a  steady  state  is  achieved.  The  presence  of  the  energetic 
plasticizer  TMETN  in  the  formulated  propellant  appears  to  increase  the 
production  of  gaseous  oxidizing  species  by  the  burning  surface.  This  in  turn 
accelerates  the  rate  of  the  exothermic  oxidation  reactions  in  the  luminous  flame 
such  that,  as  shown  by  the  CW5  propellants,  the  flame  front  moves  essentially 
to  the  burning  surface. 

The  structure  of  the  CW5  propellant  flame  suggests  that  the  oxidation 
reactions  in  the  luminous  flame  are  consuming  the  gaseous  reactants  as  soon 
as  they  are  produced  by  the  burning  surface,  thus  implying  that  the  rate  of 
condensed-phase  decomposition  controls  global  propellant  burn  rate. 

The  extinguished  burning  surfaces  of  HMX  and  RDX  CW5  and  RRHR-18 
propellants  are  shown  in  Figures  19-34.  The  extinguished  surfaces  for  the 
HMX  based  propellants  (Figs.  19-26)  clearly  indicate  a  heterogeneous  burning 
surface.  The  craters,  which  cover  the  extinguished  surface,  average  30  pm  in 
diameter,  and  match  closely  with  the  weight  median  diameter  of  26  pm  for  the 
commercial  grade  HMX  used  in  these  propellants.  As  has  been  previously 
suggested,67  our  data  indicate  that  the  burning  surface  consists  of  solid 
particles  of  HMX  resting  in  a  binder  melt  layer.  When  this  surface  was 
extinguished  via  rapid  depressurization,  the  melt  layer  solidified  and  the  HMX 
particles  burned  away  leaving  the  observed  craters.  This  observation  is  in 
contrast  with  the  extinguished  burning  surface  of  pressed  pellets  of  HMX  which 
exhibit  a  molten  layer  ranging  in  thickness  from  5  pm  at  1000  psi  to  20  pm  at 
250  psi.87 
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Although  they  exhibit  different  flame  structures  there  appears  to  be  very 
tittle  difference  between  the  burning  surface  of  the  HMX-CW5  and  HMX-RRHR- 
18  propellants.  The  only  observable  difference,  albeit  small,  is  a  slightly  more 
crystalline  solidified  binder  melt  layer  in  the  RRHR-1 8  propellant  (Fig.  26)  than 
in  the  CW5  propellant  (Fig.  22). 

The  proposed  structures  of  the  flame  and  burning  surface  for  the  HMX 
CW5  and  RRHR-1 8  propellants  are  shown  in  Figures  35a  and  b,  respectively. 
Previous  microcinematography  investigations81  of  the  burning  surface  of  HMX 
based  propellants  similar  to  our  CW5  and  RRHR-1 8  formulations  have  shown 
the  HMX  particles  to  be  coated  with  what  appears  to  be  a  molten  layer.  The 
relatively  poor  resolution  of  the  high  speed  video  system  used  in  our 
experiments  did  not  allow  us  to  confirm  this  observation.  The  results  for  the 
extinguished  HMX  propellants  indicate  that  the  HMX  remains  essentially  intact 
during  the  condensed-phase  decomposition  reactions  and  only  in  the  gas 
phase  do  its  decomposition  products  react  with  products  from  the 
decomposition  of  the  binder.  Because  of  the  apparent  melt  layer  on  the  surface 
of  the  »:rv*X  particles  it  Is  likely  that  decomposition  reactions  are  occurring  both 
in  the  solid  and  liquid  phase  during  propellant  combustion. 

The  quenched  burning  surfaces  of  the  RDX  based  propellants  (Figs.  27- 
34)  show  relatively  homogeneous  surfaces  consisting  of  a  crystalline  solid 
covered  by  varying  amounts  of  carbonaceous  residue.  The  RDX,  which 
liquefies  at  473  K  compared  to  553  K  for  HMX,  apparently  melts  aiong  with  the 
binder,  giving  a  uniform  liquid  layer.  This  is  similar  to  the  melt  layer  observed 
on  the  extinguished  surfaces  of  RDX  pressed  pellets.87  When  this  liquid 
burning  surface  is  quenched,  it  produces  the  observed  crystalline  surface. 

As  with  the  HMX,  the  RDX  CW5  and  RRHR-1 8  extinguished  surfaces  are 
essentially  the  same.  The  only  noticeable  difference  is  the  presence  of  disk 
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shaped  carbon  residue  on  the  surface  of  the  RRHR-18  propellant  (Fig.  23).  The 
amount  and  character  of  the  carbonaceous  residue,  however,  were  not 
consistent  from  sample  to  sample  and  apparently  it  is  significantly  affected  by 
the  quenching  process. 

The  proposed  surface  and  flame  structures  for  the  RDX-CW5  and  RDX- 
RRHR-18  propellants  are  shown  in  Figures  25c  and  d,  respectively.  It  !s  clear 
from  the  data  presented  that  the  decomposition  of  RDX  occurs  in  the  liquid- 
phase.  Interestingly,  even  though  the  molten  RDX  was  in  intimate  contact  with 
the  liquid  binder  there  was  no  apparent  effect  on  the  rate  controlling  processes 
as  shown  by  our  isotope  effect  studies. 

Extinguished  PB  propellant  samples  were  studied  with  an  electron 
microscope.  However,  because  of  the  extensive  amount  of  residue  on  the 
surface  of  the  propellant,  it  was  impossible  to  distinguish  between  the  burning 
surface  and  the  carbon  residue.  Previous  studies  with  an  HMX-PB 
formulation67  indicate  a  heterogeneous  burning  surface  composed  of  HMX 
particles  covered  by  a  molten  layer  sitting  in  a  binder  melt  layer.  The  data  given 
by  Beckstead  and  McCarty67  and  obtained  by  ourselves  unfortunately  do  not 
allow  us  to  determine  in  what  phase  the  HMX  decomposition  primarily  takes 
place.  Unfortunatsly,  there  is  no  information  indicating  the  character  of  the 
RDX-PB  propellant  burning  surface  as  was  the  case  with  HMX-PB.  However, 
one  can  speculate  that,  considering  the  similar  isotope  effects  with  the  CW5 
propellant,  the  surface  is  likely  to  be  a  molten  layer  of  RDX  and  binder. 

The  flame  of  the  PB  propellants  was  difficult  to  observe  because  of  the 
presence  of  sheets  of  binder  residue  which  stacked  one  on  top  of  the  other  as 
shown  in  Figure  35e.  These  wafers  of  residue  forced  the  flame  to  the  sides  of 
the  propellant  and  thus  made  it  impossible  to  determine  the  presence  or 
absence  of  a  two  stage  flame. 
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FIGURE  19.  Optical  Mircograph  of  an  extinguished 
HMX-CW5  propellant.  Magnification  is  10  times. 


FIGURE  20.  An  Electron  Micrograph  of  the  edge  of  the 
extinguished  burning  surface  from  a  HMX-CW5  propellant. 
Magnification  is  147  times.  Burning  surface  is  on 
the  left  and  unburned  propellant  is  on  the  right. 


FIGURE  21 .  An  Electron  Mircograph  of  the 
extinguished  burning  surface  from  a  HMX-CW5 
propellant.  Magni'icai'on  is  74 


FIGURE  22.  An  ."  „•  of  the 

extinguisher)  ;a.;r  ,  i ; 

Magnification 
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FIGURE  23.  Optical  Mircograph  of  an  extinguished 
HMX-RRHR-1 8  propellant.  Magnification  is  10  times. 


FIGURE  24.  An  Electron  Micrograph  of  the  edge  of 
the  extinguished  burning  surface  from  a  HMX-RRHR-18 
propellant.  Magnification  is  73  times.  Burning  surface 
is  on  the  left  and  unburned  propellant  is  on  the  right. 


FIGURE  25.  An  Electron  Mircograph  of  the 
extinguished  burning  surface  from  a  HMX-RRHR-18 
propellant.  Magnification  is  124  times. 


FIGURE  26.  An  Electron  Mircrograph  of  the 
extinguished  burning  surface  from  a  HMX-RRHR-18 
propellant.  Magnification  is  250  times. 


FIGURE  28.  An  Electron  Micrograph  of  the  edge  of 
the  extinguished  burning  surface  from  a  RDX-CW5 
propellant.  Magnification  is  121  times.  B  ming  surface 
is  on  the  top  and  unburned  propellant  is  on  the  bottom. 
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FIGURE  30  An  Electron  Mircrograph  of  the 
extinguished  burning  surface  from  Figure  17 
Magnification  is  380  times. 
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FIGURE  31.  Optical  Mircograph  ol  an  extinguished 
RDX-RRHR-1 8  propellant.  Magnification  is  10  times. 


FIGURE  32.  An  Electron  Micrograph  of  the  edge  of  the 
extinguished  burning  surface  from  a  RDX-RRHR-1 8 
propellant.  Magnification  is  122  times.  Burring  surface 
is  on  upper  left  and  unburned  propellant  is  on  lower  right. 


FIGURE  33.  An  Electron  Mircograph  of  the 
extinguished  burning  surface  from  a  RDX-RRHR-18 
propellant.  Magnification  is  136  times. 


FIGURE  34.  An  Electron  Mircrograph  of  the 
extinguished  burning  surface  from  a  RDX-RRHR-18 
propellant.  Magnification  is  500  times. 
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P  -♦  8  PHASE  TRANSITION 

HMX  in  the  solid  is  held  together  by  strong  intermolecglar  coulombic 
attractions  which  account  tor  its  high  density  and  high  melting  point.88  HMX  at 
room  temperature  is  in  the  p  crystalline  form  which  is  essentially  a  chair 
configuration  consisting  of  the  four  carbons  and  two  opposite  ring  nitrogens  in  a 
plane  with  the  remaining  nitramine  units  above  and  below  the  plane.  Above 
436  K  at  atmospheric  pressure  the  p-HMX  converts  to  the  high  temperature 
form,  5-HMX.89  The  configuration  of  8-HMX  is  best  described  as  chair-chair. 
Work  by  Brill  and  Karpowicz90  has  shown  that  the  activation  energy  and 
Arrhenius  preexponential  factor  for  the  HMX  p  -*  8  phase  transition  and  those 
for  the  thermal  decomposition  of  HMX  in  the  condensed  phase  are  essentially 
the  same  within  experimental  error.  From  their  results  they  concluded  that  the 
chemical  decomposition  of  the  HMX  appears  to  be  controlled  by  the  breakdown 
of  intermolecular  forces  in  the  condensed  phase. 

The  results  of  Brill  and  Karpowicz90  imply  that  the  KDIE  observed  by 
ourselves  and  others16-19-20  are  not  the  consequence  of  a  covalent  bond 
rupture  as  had  originally  been  concluded.  If  the  breakup  of  intermolecular 
forces  in  the  crystal  were  the  cause  of  the  KDIE  in  HMX  then  one  would  expect 
to  observe  a  significant  isotope  effect  in  the  HMX  p  -*  8  phase  transition. 
However,  this  appears  not  to  be  the  case. 

The  p  -*  8  phase  transition  for  HMX  and  HMX-d8  was  investigated  using 
differential  scanning  calorimetry  (DSC).  The  onset  of  the  HMX  and  HMX-d8 
phase  transition  endotherm  at  a  scanning  rate  of  5.00  deg/min  occurred  at 
448.9  i  1.5  K  and  447.0  ±  1.6  K  respectively.  In  addition,  the  enthalpy  of 
transition  for  HMX  was  approximately  7.5  ±  1.2  cal/g  and  that  for  HMX-d8  was 

approximately  7.0  ±  0.8  cal/g. 

The  DSC  data  qualitatively  demonstrate  the  absence  of  a  significant 


positive  isotope  effect  in  the  P  -»  5  phase  transition.  In  fact,  the  data  imply  that 
the  rate  of  the  HMX-d8  phase  transition  may  be  slightly  greater  than  that  of  HMX. 

This  observation  may  provide  some  explanation  for  the  KDIE  very  close  to  unity 
observed  for  the  decomposition  of  HMX  in  the  mixed  melt.19  It  must  be  noted 
that  the  DSC  results  represent  relatively  few  measurements  and  further  work  is 
needed  to  establish  the  actual  effect  of  deuterium  substitution  on  the  (3  -»  5 
phase  transition. 

Additionally,  if  the  conclusions  of  Brill  and  Karpowicz90  were  correct  then 
one  would  e>pect  significant  differences  in  the  intermolecular  attractions  in 
HMX  versus  HMX*d8.  This,  in  turn,  would  manifest  itself  as  differences  in  the 
crystal  structures  of  HMX  and  HMX-d8.  Recently  a  comparison  study  of  the 
crystal  structures  of  (J-HMX  and  (J-HMX-d8  was  conducted.91  It  was  found  that 
there  were  only  a  few  very  minor  differences  in  the  inter-  and  intramolecular 
distances.  These  differences  could  account  for  the  small  effect  of  deuterium 
substitution  on  the  phase  transition,  but  they  were  not  great  enough  to  generate 
isotope  effects  on  the  order  of  those  observed  in  HMX  decomposition  and 
combustion. 

From  the  data  presented,  it  would  appear  that  the  breakup  of  the 
intermolecular  forces  in  HMX  is  not  the  rate  controlling  process  in 
decomposition  and  combustion.  Furthermore,  RDX  decomposes  primarily  in 
the  liquid  phase,  and  consequently  the  breakup  of  intermolecular  forces  cannot 
be  the  rate  controlling  process.  Yet  RDX  exhibits  KDIE  of  similar  magnitude  to 
those  observed  for  the  decomposition  of  HMX.  Thus,  it  seems  unlikely  with 
compounds  as  similar  as  HMX  and  RDX  that  the  decomposition  for  one  would 
be  controlled  by  the  breakup  of  intermolecular  forces,  and  the  other  would  be 
controlled  by  a  covalent  bond  rupture. 
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CONCLUSIONS 


Our  research  has  confirmed  the  presence  of  significant  isotope  effects  in 
the  combustion  of  HMX  and  RDX.  More  importantly,  we  have  demonstrated  the 
existence  of  significant  KOIEs  in  the  combustion  of  model  nitramine  propellants. 
The  consistency  between  the  effects  of  the  pure  nitramines  and  those  for  the 
HMX  and  RDX  CW5  and  RDX-PB  propellants  indicates  similar  rate  controlling 
processes.  Assuming  the  rate  of  condensed-phase  decomposition  controls  the 
rate  of  combustion,  then  the  KDIEs  observed  for  the  pure  nitramines  and  the 
above  propellants  appear  to  be  the  result  of  a  rate  controlling  C-H  bond  rupture. 
The  combustion  of  the  HMX-PB  propellant  shows  an  isotope  effect  value 
significantly  less  than  the  pure  nitramine,  indicating  a  considerable  change  in 
the  rate  controlling  process. 

This  report  is  particularly  significant  in  extending  previous  findings  from 
pure  oxidizers  to  propellant  systems  which  are  more  realistic  with  respect  to 
additional  ingredients  which  would  likely  complicate  the  combustion 
mechanism.  In  spite  of  these  ingredients'  participation  in  the  global 
mechanism,  evidence  has  been  given  that  the  combustion  of  the  nitramine 
ingredient  in  most  cases  controls  the  global  burn  rate  of  the  model  propellant. 

These  results  represent  the  first  in-situ  studies  of  the  rate  limiting 
processes  in  nitramine  propellant  combustion.  Furthermore,  these  studies 
represent  some  of  the  first  attempts  at  applying  the  technique  of  KDIEs  to  high 
temperature  combustion.  Because  of  the  limited  experimental  data,  the 
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observed  KDIEs  can  only  be  evaluated  if  certain  fundamental  assumptions  are 
made.  Consequently,  the  assignments  of  these  KDIE  must  be  treated  with  care 
until  they  can  be  verified  by  other  methods. 
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